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EVALUATION OF DEFINITE INTEGRALS 
BY SYMBOLIC MANIPULATION* 


Abstract 


A heuristic computer program for the evaluation of real 
definite integrals of elementary functions is described. This 
program, called WANDERER (WANg's DEfinite integRal EvaluatoR) , 
evaluates many proper and improper integrals. The improper 
integrals may have a finite or infinite range of integration. 
Evaluation by contour integration and residue theory is among 
the methods used. A program called DELIMITER (DEfinitive LIMIT 
EvaluatoR) is used for the limit computations needed in evaluat- 
ing some definite integrals. DELIMITER is a heuristic program 
written for computing limits of real or complex analytic func- 
tions. For real functions of a real variable, one-sided as well 
as two-sided limits can be computed. WANDERER and DELIMITER have 
been implemented in the MACSYMA system, a symbolic and algebraic 
manipulation system being developed at Project MAC, MIT. A typi- 
cal problem in applied mathematics, namely asymptotic analysis of 
a definite integral, is solved using MACSYMA to demonstrate the 
usefulness of such a system and the facilities provided by 
WANDERER. 


*This report reproduces a thesis of the same title submitted to 
the Department of Mathematics, Massachusetts Institute of Tech~ 
nology, in partial fulfillment of the requirements for the degree 
of Doctor of Philosophy, August 1971. 
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Integrals. 


The limit concept Is bastc to mathematical analysIs. 
Being able to compute limits automatically greatly Increases 
the potential of a symbol mantpulatlton system In dotng 
analytical mathematics. In fact the evaluation of definite 
Integrals Is heavlly dependent on the lfmit process, as fs 
the expansion of functlons In power serles and many other 
mathematical problems. Although the computation of limits 
has been studied previously to some extent [11, 13], we 
describe a limftt program called DELIMITER (DEfInttive LIMIT 
EvaluatoR) whIch Is more powerful than prevlous programs. 


It Is discussed tn detall [tn chapter 2. 


The problem of computing Indefinite Integrals 
symbolically by computers has been Investigated rather 
thoroughly. First among the computer programs developed for 
this purpose was SAINT (Symbolic Automattic INTegrator) by 
Slagle tn 1961. A more powerful program named SIN 
(Symbolic !Ntegrator) [28] was developed In 1967 hy Moses. 
Theoretical work In thls area Include Richardson's 
undecidabIility result for a certaln classes of Integrals 
[22] and RIisch's dectston procedure for determining the 
existence of the Indefinite Integral of a member [n the 


class of elementary functfons (23, 24h, 25, 26]. There Is a 
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comprehensive revfew by Moses on the progress durtng the 
past ten years tn this partfcular area of symbol manlIpula~ 


tlon [21] e 


However, hardly any work has been done [n the dlrecttfon 
of definite Integration. Reported here ts a first definite 
Integratlon program called WANDERER (WANg's DEftnite 
IntegRal EvaluatoR), Experiments with DELIMITER and 
WANDERER are the princtpal subjects of this thesis. 

WANDERER has been Implemented In MACSYMA and makes use of 
many factlittes provided tn It. Some of these facilities 
are: Input and two-dimenstonal output, simplification, 
solution of polynomtal and systems of linear equations, 
canontcal rattonal functton stmplificatton and the SIN [20] 


Integration program. 


Evaluation of definite Integrals can sometimes be as 
easy as computing the Indefinite Integral of the given 
Integrand then substituting In the lImits of Itntegratton. 
WANDERER computes Integrals of thts ktnd by ustng those 
parts of SIN that have belIng Implemented In MACSYMA. 
However, many Interest!Ing definite {tntegrals are not 
obtainable [In thls manner. In some cases they may be 
Improper integrals or thelr Indefinite Integrals do not 


exist. In other cases It ts easter to evaluate the definite 


ll 


Integral directly than to obtatn the Indefinite [fntegral 
first. 


A few ltmlt and Integration problems are listed below 
to glve an Indication as to the scope of WANDERER and 
DELIMITER. The results obtaltned by the programs can be 


found tn chapters 2, 4 and 5. 


2 oo 6 
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Integrals with an Inffintte range are discussed In 
chapter 4, those with a finite range tn chapter 5. 
Important rout!Ines and algorithms are collected [tn chapter 


6. 


An application of the MACSYMA system to an asymptotic 
analysts problem Is [Illustrated [In chapter 7. The purpose 
Ts to show how such a symbol mantpulation system can be used 
to solve non-trivial problems that may occur frequently tn 
applted mathematics. This applticatlon demonstrates the 
usefulness of the many fac!ilittes provided by MACSYMA and 
WANDERER. 


A timing expertIment has been conducted to check the 
performance of WANDERER and DELIMITER, The results are 
Included In appendix C. Because of the limited character 
set of most computer consoles, some special symbols are 
needed to denote the frequently used mathematical functlons, 
constants and operators. For example, * and ** are used to 
denote multiplication and exponenttatton respectively. 


Appendix D contalns a list of notatfons we shall use, 


From Richardson's undecIdabtitty results [22], we have 


shown that the convergence of a class of Integrals of 
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elementary functlons ts recurstvely undecldable. The proof 


Is In appendtx E. 


The remainder of this chapter serves as an [Introduction 
to algebraic manIpulatton systems. If the reader ts 
famlliar with such systems he may proceed directly to the 


next chapter. 


An algebraltc manIpulation system Its a collection of 
computer programs destgned to facilitate the solutton of 
mathematical problems. Such a system has the abIlftty to 
handle both symbols and numbers. This capabllity to mant- 
pulate symbols mathematically Is what differentlates these 
systems from the various computer subroutines which 


speciallze In numertcal analysts. 


In us Ing symbol manfpulatton systems such as these, one 
usually Interacts with [ft tn a ttme-shartng environment via 
a typewrlter-IIke console. Normally, data and commands are 
typed In by the user, Results from the computer are sent 
back to the console for display. This arrangement attempts 
to provide a user the ease and flexibillty of the pencltl and 
paper he fs so accustomed to, whtle permittIng the computer 
to assist him In his algebra and formal deduction from one 


step to the next. If the manfpulations tnvolved are non- 
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trivial, so much the better. Some of the mantpulatfions such 
a system can provide are: GCD calculations, factoring, ra- 
tlonal function artthmetic, matrix mantpulatton, solution of 
algebraic equations, solution of systems of linear equa- 


tlons, tndefinitte Integratlton of elementary functtons, 


Let us take a closer look at MACSYMA as a 
representative model of other systems. MACSYMA recelves 
Inputs In the form of ltnear character strings typed by the 
user, FORTRAN=lIke notatfon Is used for the Input. For 


example, 


The user types: he (Xe*24X41)/YQ@ 


The @ slgen stgntfles the end of a command string. As a 
result of this command, a two dimenstfonal display ts 


returned. 


MACSYMA types: sm aseesecenenee 


Instde MACSYMA, expressions are represented by list 
structures In a prefix notation (common to many systems). 
For example the expression above would be represented 
Internally by somethIng ltke 

(TIMES & (PLUS 1 X (EXPT X 2)) (EXPT Y =1)) 
Inside MACSYMA. 
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Other forms of Internal representations are posslIble. 
The choice of Internal representation Is an Important aspect 
of the design of a system. Of course more than one form of 
Internal representation can be employed [tn a single 
algebraic manipulation system. In fact MACSYMA has a 
special Internal representation for ratfonal functfons whitch 
Is used to gatn efflctency during certain polynomial and ra- 


tlonal functton mantpulattons. 


One simple application of MACSYMA ts factoring the 

polynomial 

3 2 

P(X) =X + & X -11 X = 30 
A user who wants to factor P(X) using MACSYMA types: 
FACTOR(X¥¥34heXee2-11%"X-33)@ 

This Input command causes the factorizatton of P(X) over the 
Integers. The output Is 

(X + 2) (X = 3) (X + 5) 

While this problem may seem easy, factorization of 
polynomials of higher degree can be very diffftcult to do by 
hand. Indeed, algebraic manipulation systems can be most 
helpful when one wants to manipulate complicated functions 
and expressions. MACSYMA can carry out accurately with 


great speed: summation over Indices, expansion of products 
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and powers, calculation of large determinants, Inversion of 
large matrices, 

Applications [tn pure mathematics Include computattons 
In: number theory, group theory, Lle algebras and set 


theory. 


Research In the fleld of symbol and algebrafc man!Ipula-~ 
tion has ted to many new results. The fast GCD algorithm 
[4] and the finite fleld arithmetic polynomtal factoring 
algorithm [2] are two examples. ContInued work In thts area 
wlll, hopefully, result In computer systems whtch are 


Increasingly valuable to engineers and mathematIiclans. 
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CHAPTER It! 
DELIMITER 


0. Introduction 


The limit concept Is fundamental to mathematical analy- 
sis. Bas!c concepts such as the rules of differentlatton 
are dertved from limiting processes. More complex problems 
such as Improper Integratlon, convergence of serles, serles 
expansion of functtons and contour Integrals, to name a few, 
also require the computation of limits In thelr solution 


process, 


Therefore, one can expect that automatic computation of 
ltmits would greatly Increase the capability of a symbolic 
mathematics system In dotng analytical mathematics. The 
programs described below provide such a capability In the 


MACSYMA system [16]. 


Automating the computation of limits has been studied 
previously to some extent. Fenltchel [11] dlscussed certatn 
decIdabllfity problems of limits and provided, In the FAMOUS 


system, some baste routines for computing two-sided limits. 
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In his thesis, !tturrtaga [13] worked on both one and two- 
stded iimlts. In addttton to L'HospItal rule, he applied 
asymptotic analysts to quottents of polynomlals essentlally 
by replacing those polynomlfals by thelr leading terms. This 
can be done when the variable In the limit approaches Infl- 
nity so that the leading term Is asymptotically equivalent 
to the orlginal polynomltal. Moreover, he discussed limits 


of sequences. 


DELIMITER Is more powerful than elfther Fentchel's or 
Iturrtaga's limft program. L'Hosp!ital's rule Is a basic 
method used for [ndetermInate forms. In add!itlon, this pro- 
gram employs a fast routine for lfmits of ratfonal func- 
tions. It has an effictent algorithm for a class of expres- 
stons called RP~expresstions which tnvolve radicals of poly- 
nomials. It also applfes the method of reductng complicated 
expresstons by replactIng subexpresstons with asymptotically 
equivalent expressions. The method of compartng orders of 
Infintty of expresstons and several other heurfstIic methods 
are used. In some cases, power serles expansions are employ- 


ed to obtain the lfmits. 


This chapter [fs based on a paper by the author [29] 
presented at The Second Symposfum on Symbolic and Algebrale 


Mantpulatlon, Los Angeles, March, 1971. 
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1. Deftnittons and Symbols 


DELIMITER ts destgned for finding the limit of a 
single-valued function f(X) of a real or complex varlable X 
as X approaches some limit potnt. One-stded limfts can be 
computed If X and f(X) are real-valued. The classes of 
functions allowed by this program Include ratlonal, radical, 
logartthmic and exponential functtons, and also the trigono- 
metric functtons SIN, COS and TAN, and the hyperboltc 
functions SINH, COSH and TANH. The general form of a user 


command Is 


LIMITCexpresston,vartable,value,directton)@ 


with the fourth argument optlonal. PLUS as a fourth argu- 
ment indIcates the one-sided limit from above, MINUS from 
below. The absence of the fourth argument Indicates no 
restriction on the dtrectton of approach. One~stded limits 
are not allowed for complex-valued functions or limft 


pofnts. 


If F(X) Is not continuous at X=a, the two-sided limit, 
LIMITCF(CX),X,a), does not exist. In such a case the symbol 


UND fs returned by the program as an answer. 
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The following are some symbols used by this program. 


ZERO+ 0+ ZERO= 0- 

UND undef Ined 

IND Indefinite but bounded 

INFINITY the poInt at Infinity tn the complex plane 


INF and MINF whll conttnue to Indicate respectively 
positive and negative Infinity. Here are some examples of 


how these symbols are used Internally. 
LIMITCCOS(X),X,%P!1) = -1+ ZERO+ 
LIMITCSINCX),X, INF) = END 
LIMITC1/SINCX),X,INF) = UND 


LIMETCI/(X=21),X,31) = INFINITY 
The direction from which a limit pofnt ts approached fs 
Important. By use of ZERO+ and ZERO-, simpltfication rules 
such as 
1/(-1 + ZERO+ + 1) ——»> INE 


and 
1/(-1 + ZERO- + 1) —» MINE 


are possible In DELIMITER, 
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2. Baste Rules and Outline of Algortthm 


Some of the rules used In thls program are the so-call- 
ed "trivial" ones for limits of continuous functtons, 
namely, the limit of a sum Is the sum of the limits, etc. 
There ts a limft routtne for each of the functtons SIN, COS, 
TAN, SINH, COSH, TANH and LOG (base %E). For other func 


tions the rule 


LIMITCFCe(X)),X,L) = LIMITCFCY),Y,LIMITC g(X),X,L)) 


Is used. 
For one-sided limits where the varltable approaches a 
potnt other than 0, a change of variable ts made to bring 


the polnt to the orfgtn. For example, the limftt 


LIMITCF(X),X,a, PLUS) 
Is converted to 


LIMITCFCatY),Y,ZERO+) 


There Is, of course, a complete set of rules governing 
simplification of the new symbols. They are all of the fol- 


lowlng nature 


22 


IND/INF = 0 -1/INF = ZERO- 

~ ZERO+ = ZERO- ~ INF = MINF 
1/INFINITY = 0 “1/INFINITY = 0 
MINFe*2 = INF (ZERO+)+ZERO=- = 0 


LOG(C+ZERO+) = LOG(C) + ZERO+ 
We wlll now briefly outline the algorithm, and then go 
Into the detalls of some of the component routtnes. Let us 


consider 


LIMIT CEC(X),X,L) 


where L Is any number or symbol Including INFINITY, INF, 
ZERO+ and ZERO-. If b= MINF, [t ts set to INF by a change 
of variable Y=-X. Upon recefving the arguments, namely 
E(X), X and L, DELIMITER checks whether L=INFINITY or E(X) 
Involves “I. If elther or both Is true then a global 
Indicator CPLX Is set to the value TRUE whIch [Indicates that 
the given limit Is to be evaluated over the fleld of complex 
numbers. Otherwise CPLX Is set to the value FALSE which 
slgntfles a limit problem over the reals. If CPLX=TRUE, all 
use of INF and MINF are replaced by INFINITY and the notion 
of approaching a limft potnt from one side Is no longer 


valld. 
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The followtng diagram Its a stmplifted overview of 


DELIMITER. 
DELIMITER 
SIMPLE LIMI T————--——>RETURN ANSWER 
RATIONAL FUNCTION=——=>RETURN ANSWER 
Rie CO pasar RE IURE ANSWER 
BASIC LIMIT L'HOSPITAL'S SERIES ORDER 
RULES RULE EXPANSION OF INFINITY 


Fig, 2 


As indtcated In fig. 1, DELIMITER has four stages. 
Results of simple limit problems are Itmmedtately returned by 
the first stage. Simple limits Include two cases: 1) E Is 
Independent of X. 2) E(X) = X. In the second stage, a fast 
algorithm Is used to obtatn limits of ratfonal functions. 
The third stage employs an algorithm for limits of RP-forms. 
If the first three stages do not produce an answer, the 
given limit problem enters the fourth stage which contains a 


varlety of methods Including the four principal ones shown 
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In the figure. In thls stage the program willl first try to 
extract the numerator N(X) and the denominator D(X) of E(X) 
so that E(X) = N(X)/D(X). If both NCX) and D(X) are RP= 
forms the limit Is computed by an algorithm spectally for 
these forms (see section 5 for definition and algorithm). 


Otherwise the rule 


Tim ECX) = lim NOX) /1tm DCX) 


Is applfled. If both lim N(X) and Im D(X) are 0 or Infinite 
we have an Indeterminate form. In this case L'Hospital's 
rule will be applted to E(X) with one exception: when L fs 
INF and lim N(X) and lfm D(X) are both Infinite and both 
N(X) and D(X) contain exponential functions of X which tend 
to INF as X approaches L, In this case the method of com= 
paring orders of InfinIity, to be descrftbed shortly, fs ap- 
plied. If DELIMITER can not find an answer Indication of 


fallure wlll be returned. 


Sometimes the program needs to know the value range of 
a symbolfc parameter In order to compute the limit. In such 


cases, the program wlll query the user at his console. 


It Is Important to note that these programs do not 
store any table of limits Thus every lImfit obtalined ts a 


result of computatton. 
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Other methods available [Include : change of vartable, 
discontinulty tests and analysis of the behavior of a func- 
tion near a finite point. <A flowchart [tn append!Ix A pre- 


sents the flow of control In a more complete manner. 
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If n=0, the answer Is c;} 


it) 


tf n < 0, then Tf CLPX TRUE, answer fs 0, but If 
CPLX ts FALSE then the answer —[s ZERO+ provided 
the stgn of cxlLe*n Ts + and ZERO- [tn case the slgn 
of c*elL*een ts =3 

otherwise, (n > 0) the answer Is INFINITY tf 
CPLX=TRUE, and the answer Is INF or MINF, 


depending on the stgn of c*lL*e*n, Tf CPLX=FALSE, 


28 
4h. The Limit of a Quotlent of RP-forms 


An RP-form Is deffned as an expression obtainable by 
combInIng polynomtals and posttIive rational number powers of 
polynomials ustng the operators +, -, *, Here are some 
examples of RP-forms, 

3 1/3 
F(X) =SQRT(C4 X +5 X) + (X +1) +X +C 
2/3 

e(X) = AX SQRTC(X+1) + X 

Note that the definition of RP-forms does not allow 
nested radicals of X. For Instance the follow!Ing Is not an 


RP=-form. 


2/3 
(X + SQRT(X)) 


Now constder the limit problem, 


LIMITCE(X),X,L), ECX) =N(X)/D(X) 


where N(X) and D(X)#0 are RP-forms and at least one of them 


ts Irrational. Note that N(X) or D(X) can be a constant. 


Let us define the operators EXPO and COEF by 
EXPO(N(X) ) 

=hfighest exponent of X tn N(X), 
COEFCE,N(X)) 


=the coefflctent of E In N(X). 
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For example, 


EXPO(g(X)) =7/6, 
COEF (X**(3/2),f(X)) =2 


EXPO fs obtained by a sImple special purpose routine; COEF 
1s computed by ustng the RATCOEF routine In the ratfonal 


functton package of MACSYMA [16]. 


Now the limlt {ts computed by the followIng algortthm. 


(1) If L =MINF, X Is replaced by -X, L by INF. 


(2) (If L =INF, the follow!Ing asymptotic analysIs ts 


made. 


1) compute a =EXPO(N(X)), b=EXPO(D(X)). 


Tt) Let N1(X) and D1(X) be NCX) and D(X) with 
polynomftals under radicals replaced by thelr 


leading terms, respectively. Compute 


a'=EXPO(N1), b'=EXPO(D1). 


Tit) If efther a#a' or b#b', then ft ts usually 


suffictent to replace each radical In E(X) by 
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D(X) before ustng L'Hospital's Rule. 


If NCL)=0 and D(L)#0, 
the answer Is 0 If CPLX=TRUE, otherwise the be- 
havior of D(L)*#N(X) at L ts examtIned to decide 


whether 0, ZERO+ or ZERO- {s the answer. 


If D(L)=0 and NCL)#0, 

the answer Is INFINITY If CPLX =TRUE, otherwise 
the behavior of N(L)*D(X) at L Is examined to 
decide whether INF or MINF [fs the answer. 


Otherwise, the answer Its ECL). 


For large X the serfes expanston 


2 1 
SQRTCX + 1) = SQRT(-- + 1) X 
2 
X 
1 1 
= X (1 tS ey oe SSS + . .) 
Z 4 
2 X 8 X 
1 1 
= Xx + —— = asa cali + ee e e 
2 X 3 
8 X 


is convergent. It Is clear that this expansfon method can 
be applled to an arbltrary polynomial to a positive frac~ 
tlonal power. This fact establtshes the validity of step 


(2-Ttl). 
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5. Indeterminate Forms and L'Hospital's Rule 


L'HospIital's rule {ts applfed to an Indeterminate form 
E(X) of the form 0/0 or INFINITY/INFINITY. Ut flrst 
calculates N'(X) and D'(X), then, 


If D'(X)=0, L'Hospital's rule falls, otherwlse, simplify 


N'(X)/D'(X) and evaluate 


LIMIT (N'CX)/D'(X),X,L). 


The expresston N'(X)/D'(X) can be more complex than 
E(X) and successfve application of this rule may lead 
nowhere. Therefore the number of times this rule ts called 
successIively fs counted and the relative complexity of 
N'(X)/D'(X) to EC(X) Its tested to decide whether to continue 
this approach. Our criterfon of complexity Is based on the 
number of distinct nonratfonal components of an expression. 
If this number grows for three consecutive times, the ap- 


plicatton of L'Hospital's rule fs halted. 


The tndeterminate form O*INFINITY Its transformed to 
either O/C1/INFINITY) or INFINITY/(1/0), depending on which 
Is simpler, before applytng L'Hospital's rule. Other In- 
determinate forms such as le*INFINITY, INF**0 and O**0 are 


handled by the logarithmic reductton : 
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F(X) **G(X) = ZEx*( G(X) *LOGCF(X))). 


Theoretically the [tndeterminate form (INF = INF) can 


always be reduced to 0/0 by rewrltting [ft as 
(1) C1/INF = 1/INF)/CI1/INFeINE), 


but thls method often makes the expressfon much more 
complicated. It Its useful, though, for expresstons 
Involving trigonometric functions as can be seen tn 


1 1 
LIMITC= = ------,X,0,PLUS) = 0. 
X SINCX) 


DELIMITER transforms the given expression In this example to 


SINCX) =X 

X SINCX), 
then appltes L'Hospital's Rule to obtain the answer 0. 
Therefore, the method (1) Is used for expresstons Involving 
trigonometric functions. Other types of expressions can be 


dealt with more readily by comparing degrees of Infinity of 


the subexpresstons or by serles expanston. 
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6. On the Sten of Infinity 


If a functlon f(X) Is discontinuous at X=a, another 
function g(f(X)) Is not necessarily discontfnuous at a. 
For example 


1/x A 
LIMITCCA X + 1) X%,0)2 ZE 


Therefore, discontinulty may be encountered [tn the course 
of computing a limit of a continuous function. Constder 


I/ECX) with ECa)=0 and the limit problem 


LIMITC1/E(X),X,a), 


lf CPLX = TRUE, the answer [ts INFINITY. Otherwlse, a rou- 
tine named BEHAVIOR is used to analyze the behavior of E(X) 
near X=a. If E(X) approaches 0 from above or below as X ap- 
proaches a, then the answer Is INF or MINF, respectively. 


Otherwise, the answer [ts INFINITY. 


The BEHAVIOR routine uses differentiation to analyze 
the behavior of a function near a pofnt by Investigating Its 
slope or concavity at the pofnt. This routIne Is also used 
by other programs. One of these ts the program whtch com- 
putes limits of the trigonometric function tangent. In case 


the argument of TAN approaches %PI/2, say, tt fis Important 
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to know the direction of approach. 


7. Heurtsttcs for Compartng Orders of Infinity 


Let f(X) and g(X) be two real-valued functions which 
become positively Infinite as X approaches INF. We define 


the symbol >> by 
Definition : g(X) >> F(X) If 
LIMITCFCX)/20X),X,INF) =0. 


Ustng thts symbol we can make the followIng brief 
table. 


wae SERRZER*XD >ZERXD> OX ¥en? LORZ(X) > 1og(log(X)) we 


Many limit problems with the varlable approaching INF 
can be solved very effictently by using thls concept. Some 
of these problems can be difficult to solve by other me- 
thods. For example, L'Hospital rule falls to compute the 


limit 
LIMITCCZE+1) eaXen2/ZEweX,X, INF), 


while the answer Its obviously INF. We Incorporate this con- 
cept by a routine, STRENGTH. It can classify the order of 


infinity of the argument according to the follow!ng rules, 
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STRENGTH(c)=(c) 
STRENGTH(LOG(R(X)))= LOG 
STRENGTH(A**B(X))= EXP 
STRENGTH( P(X) **(m/n))= med/n 


where c Is [ndependent of X, m and n are positive [ntegers, 
d the degree of the polynomial P(X). Note the STRENGTH of a 
constant Is denoted by that constant Inside parenthestIs and 
the STRENGTH of a polynomtal [fs a constant. BCX) and R(X) 
can be any functfons In X, except those which can cause the 
relevant argument of STRENGTH to be sImpliffed Into one of 
the other three cases. \le can assume A to be %E for If A Is 
any other expressfon It Its always possible to change the 
base to %£. STRENGTH of a sum Its the maximum of the 


strengths of the terms [In the sum. 


There Is a bastc compartson rout!Ine which knows the re- 


lations between the strengths of functions, 
EXP>>a>>LOG>>(b) and c>>d If cod 
where a, b, c, and d are constants. 
Let N(X) and D(X) be two products satisfying 


LIMITONCX),X, INF) =0NF 
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and 


LIMITCD(X),X, INF) SINF 


and no factor In N(X) or D(X) Is a sum or a trigonometric 
function. To compare the relative order of Infinity of 
NCX) and D(X), the followIng method ts used. The answer 1 
Indicates N(X)>>D(X), -1 Indicates D(X)>>N(X), 0 no 


dectston. 


Algorithm COMPARE(N(X),D(X)): 


1) Remove from N(X) and D(X) any common factors. 


2) Apply the STRENGTH routtne to each factor of N{X) 


and obtain the maximum strength SN. 


3) Do the same to D(X) and obtatn Its maximum strength 


SD. 


k) If SN>>SD the answer fs 1, If SD<<SN the answer [s - 
l. 


Otherwlse, 


1) Let 
N1(X)=product of all factors with strength SN [In 
N(X). 


D1(X)=product of all factors with strength SD [In 
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D(X). 


it) If SN=SD=some constant, set 


SN=STRENGTH(N1) and SD=STRENGTH(D1). 


If SN>>SD, the answer fs 1, If SN<<SD the answer Is 
“1. If both N1(X) and D1(X) are polynomials In X, 
the answer Is the leadIng coefftctent of the 
polynomtal N1(X) - D1(X). Otherwise If the varlable 
INDICATOR has the value 1, this means each of N(X) 
and D(X) Is an exponent of an exponent of an expo- 
nential expresston and they come as a result of 
recursive calls to COMPARE. For example N(X) may 
come from ZE*x*ZE**eN(X). In this case, the algorithm 
goes as follows. First INDICATOR ts set to 0 then 
E=LIMITCNCX)=-D(X),X, INF) Is computed and the answer 
Is 1 If E>0, -1 If E<0, 0 otherwise. 


If SN=SD=EXP, 


The exponents A(X) and BCX) of N1(X) and D1(X) 
are computed, That Is N1(X) =%E**A(X) and D1(X) 
=%ZE**B(X). If both A(X) and B(X) are polynomials In 
X, the answer {fs the leading coefficfent of the 
polynomial ACX) - BCX) If this polynomfal [fs not a 


constant, 0 If it Is. If A or B fs a sum, A {ts set 
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to LIMITCACX)-B(X),X,INF) and the answer Is 1 If 
A=INF, -1 If A=MINF, 0 otherwise. In case nelther A 
nor B Is a sum, If both A and B are of strength EXP, 
set the varlable INDICATOR to 1. The answer Is ob- 
tatned by evaluating COMPARE(A1(X),B1(X)), where 


A(X) = ZEe*Al(x) and B(x) = ZE**Bl(X). 


If SN=SD=LOG, 


o 


If N1(X) Is LOGCFC(X)) and D1(X) Is LOGCg(X)) then 
set SN to STRENGTH(f(X)), SD to STRENGTH(g(X)). Now 
the answer Is 0 If both SN and SD are constants, 1 


If SN>>SD, -1 If SN<<SD. The answer [ts 0 otherwise. 


We shall next dIscuss the [ItndetermInate form (INF-INF) 


and see how COMPARE can be used In such situatfons. Let 


where 


7 
F(x)= Do f (xX) 
zh 


LIMITCF, (X),X, INF) =|NF or MINF 


for I=1,2,...n. 


LIMITCFCX),X,INF) becomes Indeterminate when there exists |} 


and j such that 


and 


INF 


LIMITCF, (X),X, INF) 


LIMITCE (X),X, INF) MINF 


40 


The method of orders of Infinity can sometimes be ap- 


plied to this probiem. The algorlthm used fs as follows. 


The program forms two lists L1 and L2 of the f,'s whfch 
go te INF and MINF respectively. Then from each of L1 and 
L2 a member of maxImum order of Infinity can be obtalned, 
Let these be r(X) in L1 and s(X) In L2. Now if r(X) >>s(X) 
the answer Is INF, If s(€X)>>r(X) the answer fis MINF,. Other- 
wise, the problem can be very difficult. However the pro- 


gram trles to evaluate 
INFeLIMITCECX)/s(X),X, INF), 
which sometimes producés an answer, 


Now we will follow the major steps of the solution of a 
problem In a more detalled manner, Let 
2 


X*SQRTCX + 1) 
A(X) 


u 
oo 
m 


BCX) = ZE 
Consider the problem LIMITCACX)=-B(X),X,INF). First the 
Indeterminate form (INF-INF) Is encountered through the fol- 


lowIng steps. 


LIMITCACX)-BCX),X, INF) 


kl] 
LIMITCACX),X,INF)-LIMITCBCX) ,X, INF) 
LIMITCACX),X, INF) =%E*e*INF =INF 
LIMITCBCX),X,INF) =%E**INF =INF 


INF=INF 
Then COMPARE ts applted to A(X) and B(X), only to 
arrive at no conclustfon about the relative order of Infinity 


of A(X) and BCX). 

COMPARE (CA(X),BCX)) 
STRENGTHCACX) ) = EXP 
STRENGTH( BCX) )=EXP 

2 2 
COMPARE(CX*SQRT(X + 1),X ) 

2 
STRENGTHOCX*SQRT(CX + 1)) = 2 
2 

STRENGTH(X ) = 2 


Now the problem Its converted to the following form and the 


answer Is INF, 
INFeCLIMITCACX) /BCX),X, INF) -21) 


2 2 
LIMITCX*SQRTCX + 1) - X ,X,INF) = 1/2 
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INFe(ZEee(1/2)-1) 


INF 


8. Examples 


A number of IfmIt problems solved by DELIMITER are 
included In this section. They are presented In the form of 
actual [nputs and outputs of the MACSYMA system. Lines 
labelled (Ci) are !nput or command lines and (DI) output 


lines or answers. 


(C1) LIMITCX**LOG(1/X), X, INF)Q 
(D1) 0 
(C2) (COS(X)-1)/€ ZE*xXe*2 -1)@ 


cos(x) - 1 
(D2) 0 ee en 
2 
xX 
ZE - 1] 
(C3) LIMIT (D2,X,0)@ 
1 
(D3) was 
2 
(Ch) (C1#A*X)**(1/X)@ 
1/X 
(D4) (A X + 1) 
(C5) LIMITCD4,X,0)@ 
A 
(D5) %E 


(C6) X¥*2aClaXeeht5 ee (1/2) -24XeahG 


(C6) 


(D6) 
(C7) 
(D7) 
(C8) 
(D8) 
(C9) 
(D9) 
(C10) 


(D10) 
(C11) 
(D011) 
(C12) 


(012) 


(C13) 
(D13) 
(C14) 
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Xew2e(haXee 5) ¥# (1/2) -24XeehQ 


2 4 4 
X SQRT(K X + 5) = 2 X 


LIMIT(D6,X,INF)Q 


5 
4 
L/X-1/SINCX)@ 
1 1 
Xx SINCX) 
LIMIT(D8,X,0,PLUS)@ 
0 


Ewa (Xe (XH*24] ee (1/2) ) -ZEweXee2G 
2 
X SQRT(X + 1) X 
ZE - ZE 
LIMIT(D1,X,INF)@ 
INF 
(ZE*eX+X¥LOG(X))/CZERR CXR 3541) 4401/2) +LOGCX* eh 4X41) )@ 


3 
rf SQRT(X + 1) 
LOG(X +X + 1) + ZE 


LIMITCD4,X,INFI@ 
0 
1/(X**35-64X ee 2+114*X-6)G 


(Dib) 


(D15) 


(C16) 


(D16) 


(O27) 


(D17) 


(C18) 


(018) 


(C19) 


(D19) 


Oe ee er ee ee ey 


LIMITCDI4,X,2,M!INUS)G 


INF 


(X*SQRTCX+5) 41) /(SQRTCU*X ee 3541) +X) C 


a SQRTCX + 5) + 1 


ee ee ee ee! 


SQRTC4 X + 1) + X 


LIMITC016,X,1INF)G@ 
1 


2 


TAN(CX)/LOGCCOS(X) DG 


TANCX) 


ee ee ee ee) 


LOGCCOS(X)) 


LIMETCD1I8,k%,ZP1/2,MINUS)@ 


MINE 


bk 


(C20) 


(D020) 


ee en ee ae 


LIMCOEFF CD20, 


£),2,%1*5P1/20G 
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CHAPTER IT! 
WANDERER=AN INTRODUCTION AND OUTLINE 


1. Introductton 


The present chapter and the three subsequent ones are 
devoted to a complete description of WANDERER, a heurlstIic 


program for evaluating definite ftntegrals over a real range. 


WANDERER cannot solve all definite Integration 
problems, but [ft can certainly obtaln solutions to a large 
number of Integrals. Many examples presented are actual 
problems taken out of graduate text books [{1, 5]. The basIc 
approach, In WANDERER, for evaluating an fntegral Is by com- 
plex contour Integration and residue techniques, This ts a 
very general method particularly useful In evaluation of 
Infinite integrals. Some other methods avallable to 
WANDERER are: substitution, Integration by parts, diffe- 
rentiatlon with respect to a parameter, pattern recognition 
and table look-up, finlte-to-Infinite conversion, [ntroduc- 
tlon of a parameter and partition and transformation of the 


range of Integration. 


47 


The user's command In MACSYMA for definite Integration 


Is the key word DEFINT. DEFINT takes four arguments as [n 


DEFINTCexp, var, a, b) 


where exp Is the Integrand, var the vartable of Integration, 
a and b the lower and upper limits of [tntegratton. There Is 
no restrictlon as to what types of [tntegrands are allowed as 
input, as long as they remafn elementary functions. 

However, WANDERER requlres the ftntegrand to be finite [In the 


Integration range except possibly at the end poftnts a and b. 


A few special symbols are used In WANDERER for some 
branches of the logartthm function. They are listed In 


append!x D. 


2. Outline of Approach 


Before going Into the detalls of methods and algorithms 
for evaluation of the many types of definite Integrals whitch 
wlll be discussed In the next three chapters, a quick look 
at the whole picture with emphasIs on the flow of control Is 


In order. 


WANDERER Is a heurlstic program whitch computes definite 


Integrals by trying to apply one or more of the methods or 


algorithms built [Into It. The clues as to which of the 
methods to apply are obtained by examfintng the range of 


Integration and the form of the integrand. 


Fig, 1 on the next page serves as a simple outline of 


the flow of control in WANDERER, 
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SS) 


WANDERER 
S=success 
F=fallure 
stmple Itntegrals—ereturn answer 
F 
normalfzatton 
of tntegrand and 
Integration range 
methods for methods for 
Infinite tntegrals, finite Integrals, 
applted according to Including proper and 
range of Integratton improper Integrals 


more general methods 


return answer transformatlon, ———»return answer 
change of varlable 


4 


divergence test 


Fle, 1 


As Indteated tn flg. 1. WANDERER can be said to have 
five stages. In the first stage simple Integrals are com- 
puted whIch t{nclude cases such as a constant [Integrand, an 


Integral with equal upper and lower limits of [ntegratton 
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and, If the varlable of Integratton Its X, Integral with 
integrands which are polynomlals tn X or ZE**X. The 
normallzatton stage wlll make sure, by transformation If 
necessary, that (1) the lower limit fs less than the upper 
limit, and (2) constant factors, tf any, In the Integrand 
are removed to be mult!tplted tnto the final answer of the 


Integral. 


The third stage comprises most of the methods that wlll 
be described in chapters 4 and 5 and Is the work-horse of 
the whole program. Many methods In thIs stage are grouped 
according to the Integration range In which they are 
appropriate. Finite Integrals of ratlonal functions are 
transformed Into Infinite Integrals which will be evaluated. 
For other finite Integrals, an effort Its first made to 
determine whether they are [tmproper Integrals. The absolute 
divergence of a finite Improper tntegral Is tested before 
any attempt at evaluation. If the third stage fails to 
obtaln any results, the given problem enters the fourth 
stage which contains more general methods approprtate to 
varlous types of Integrands. A transformation or change of 
varlable Is often done In this stage. When WANDERER runs 
out of methods, the convergence of the given tntegral will 


be challenged, In case tt Is divergent, WANDERER will so 
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Indicate In the output. Otherwise, [IndIcatton of failure 
will be returned. I!t may be argued that tests for 
convergence should be conducted before evaluation, but this 
Is not necessary because WANDERER Its so destgned that 
whenever {t produces an answer, the given [Integral Is 
convergent. That Is to say each Individual method has fts 


own convergence and divergence cond!ittons bullt [n. 


The limit program described In the previous chapter Is 
used whenever a limit computation Is needed In WANDERER. It 
may often be requfred In changing the varlable of Integra- 
tlon, computing restdues and testing for convergence or 


divergence. 


3. Outline of the Computation of Restdues 


One of the most powerful methods tn evaluation of 
definite Integrals is the use of contour fntegration through 
residue calculations. The [mportance of contour [ntegration 
and the residue theory [fn the sequel warrants a brief 


summary of relevant facts from complex analysts [5]. 


(i) The Cauchy !ntegral Theorem: Let D be a simply 
connected domain and let f(Z) be analytic In D. Let C be a 


closed contour fn D. Then 
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[ea dZ = 0 
¢ 


A simply connected domatn ts Intultively an open set of 
polnts with no holes or cuts In Its Interfor. For example 
the disk [Z| < 0 Is stmply connected whIle the annulus 1 <¢ 
IZ{ < 2 Is not. The Cauchy Integral theorem ts one of the 
most Important tools In complex analysts. It Is by use of 
this theorem that path of [ntegratlon can be deformed. More 
precisely, If Cl and C2 are two different curves [n D lead- 
Ing from the same starting potnt v to the same end potfnt w, 


as shown In flg. 2, 


Fils, 2: 
and f(Z) Is analytic In D, then 


[r@ d2 = fs dZ. 
dt C2 


Therefore, the contour Cl can be deformed tnto C2. The 
above Its an Immedtate consequence of Cauchy's Integral 


theorem which gives 
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[rw dZ - fee dZ = 0 
ci C2 


(tl) Definition of Restdue: Let f(Z) be analytic In the 
punctured nelghborhood of a potnt p. Let Cp be a small 
positively orlented circle with center at p. Then the 


restdue of f at p ts defined as 


1 
si fe dZ 
221 

Cp 


The value of this Integral ts the coefficlent of the 
term 
-1 
(Z = p) 
In the Laurent expanston of f about the potnt p, so this 


coefflclent can also be regarded as the residue of f at p. 


(lll) The Restdue Theorem: Let C be a simple closed 
positively ortented contour. Let D be a simply connected 
domatn contalnitng C and its [Interfor. Let Al, A2, . « «, An 
be points Inside C. Except for Isolated stngularities at 
Al, A2, . . e, An, let f(Z) be an analytic functlton In D0. 


Then 
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n 
fa dZ = 2%t } (Restdue of f at Aj) 
c Jet 

The algorithm for restdue computatfon used tn WANDERER 
wlll be deserfbed briefly here. Detalled dIiscusstfon of thls 


algorithm can be found tn sect. 6-1. 


Let f(Z) be a functfon of a complex variable Z, 
analytic everywhere In a domain D except for a number of 
poles In D. Suppose F(Z) can be written fn the form 

F(Z) = UCZ)/V(Z) 
such that U(Z) fs analytic In D. This means, that poles of 
f(Z) are zeroes of V(Z). Suppose p Is a pole of order m of 
F(Z), the residue of f(Z) at p Is computed by the followIng 


algorithm, 
RESIDUE ALGORITHM: 


If m = 1, compute as answer 
UCp)/V'(p) 
otherwise, If V fs a polynomial, 
(1) Set V to the quotient of V(Z)/(Z - p)**m 
which Its computed by long diviston. 


(if) Return the result computed from 


m=-1 
1 d U(Z) 
(m - 1)! (3 V(Z) 


Z=p 
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otherwlse, (m > 1, V not a polynomial) compute and 
return as answer the lfmit 


m-1 


1 d m 
Lim — (Z - p) F(Z) 
Z~p} (m - 1)! \dZ 


DELIMITER ts used In obtafntng such a limit. Methods 


for infinite Integrals are dIscussed In the next chapter. 
Those for finite tntegrals are !tncluded tn chapter 5. 
Important algorithms and algortthms common to many methods 


are detailed In chapter 6. 


56 


CHAPTER IV 
SYMBOLIC COMPUTATION 
OF INFINITE INTEGRALS 


0. Introductton 


The alm of this chapter ts to give a detalled account 
of the methods employed by WANDERER for Improper Integrals 
with an Inffnite range. Contour Integration and residue 
computatton play a very Important role tn many of the 
methods. The algorithm for Integrals with rattonal 
Integrands, which wlll be discussed first, ts the most com- 
Plete. Sectlons are formed according to the function types 
of the Integrand and ordered roughly tn [tncreastng com- 
plexity. Quite a few examples are Included, some of them 


with references to books or Integral tables Indicated. 
l. tnfintte Integral of a Ratlonal Functton 


In thts sectlon methods for evaluating Infinite Inte- 
grals of a rational function wlll be presented. The range 
of Integratton Its from elther 0 or minus Infinity to 


Inffnity. Other possible Infinite ranges such as (a INF), 
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(MINF a), where a ts finftte, can be converted to (0 INF) by 
a change of varlable [In the Integral. 


First let us consider 


co 
L -[ R(X) dX, ROX) = PCX)/Q(X) 
-00 


where P(X) and Q(X) are polynomials tn X over the fleld of 
complex numbers. 

WANDERER requires deg(Q(X))-deg(P(X)) 2 2 to Insure the 
convergence of the given Integral. If R(Z) has no real 


poles, then The integral L can be computed by evaluating the 


J .[ R(Z) dZ 
C 


around a famlilfar semt-ctrcular contour tn the upper complex 


contour Integral 


Z-plane (flg. 1). One can eastly prove that L = Jas fT 


tends to Infinity. 
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Therefore 


(1) L = 2 %PI 31¢ y. Res(R) over poles of R Inside C) 


The restdue subroutine outlined In the previous chapter Is 
used to compute the resfdues needed In this formula. If 
R(Z) has poles on the real axis then L ts divergent. 
However, {tf the real poles of R are all simple then the 
Cauchy PrinclIpal Value of the Integral L exIsts, and can be 
obtalned by Indenting the contour C at these singular potnts 
on the real axis. For such Integrands WANDERER will compute 


this principal value, which ts given by 
(2) (P)L = 2 %PI %I( ) Rés(R) over poles of R Inside C) 


+ ZP1 21€ Res (R) over sImple real poles of R), 


as an answer to the Integral L. (¢(P)L stands for the 
princtpal value of L. Whenever the answer Is a principal 


value, the message PRINCIPAL will be sent to the user first. 


Here [ts an example solved by WANDERER 
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(C1) (X¥*2+A¥X4+B)/(X¥*b +10 eX 44249) Q 


(Dl) «—_ tee ee ene 


(D2) eee ee 


The expresston Dl has four simple poles: X = %l, ~4l, 
3 1 and -3 31. Ustng formula (1), only residues at X=%I 
and X=3 %1 need be computed. The residues are 
at X=l (B-1+A %1)/(16° 31) 
at X=3 Z| ~(B=9+3 A Z1)/C48 31) 


These values are computed by the residue algorithm given In 


the prevflous chapter. 


The full algortthm depends on findtng the poles of R. 


How thts fs done wlll be discussed tn detall after constder- 


o 
K + [rw dX 
°o 


if R(Z) has no real pole which Is positive or zero, 


Ing the next Integral. 


this Integral can be evaluated by Integrating the following 


contour [Integral [30] 
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J -f f(Z) dZ, f(Z) = PLOG(=Z) R(Z) 
c 


where PLOG denotes the principal branch of the LOG functton, 
i.e., the imaginary part of LOG(-Z) lfles between %PI and = 

Pt. The contour C cons{ists of circular arcs of radtl a, b, 
and two straight lf!nes jolIning thelr end-points as shown In 


fig. 2. 


One can verify that contrfibutlons from the two clircular 


arcs of C vanishes In the lfmft. Thus, 


w 6 2 SPL &t 
J -{ PLOG(-X) R(X) dX + i PLOG(-X 3E ) R(X) dX 
t-) 


@ 


o co 


(+) 
- { PLoscx) R(X) dX -{ PLogcx) R(X) dX = 2 SPI zif ROX dX 
° ° ° 
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Hence the value of the Integral K fs given by 


(3) K = “(> Res(f) over poles of R ftnside .C) 


Therefore, the value of K ts -1 times the sum of 
residues of LOG(=<Z)R(Z) In the complex Z-plane cut along the 


positive real axts. 


It may be of some Interest here to mention that a 
theorem in |30] gives a value for the tntegral K which 


differs from the correct formula (3) by a sign. 


Now suppose R(Z) has poles which are real and positive, 
then the Integral K ts divergent. But tts Cauchy principal 
value exists If every such pole Is of order 1. This value 
is given by -1 tImes the sum of residues of LOG(-Z)R(Z) In 


the entire Z-plane punctured at Z=0, as given fn 
(4) (P)K = -( )) Res(f) over poles of R). 
K diverges if R(Z) has a pole at Z=0. The following 


examples were computed by WANDERER. For ease of reference, 


expression (D1) Is repeated here 
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X + AX +B 
(D1) ee en nee nen 


k 
X +10X +9 
(C3) DEFINT(D1,X,0,INF)@ 
2 ZP! B + 6 LOG(3) A + 6 &PI 
(D3) ee meme ene enn 
The answer [tn (D3) ts obtatned by locattng the four 


stmple poles of (D1) and applytng formula (3). 


(C4) 1/(X*¥*¥24+X+1)@ 
1 

(D4) eames eencne-= 

2 

(X +X + 1) 
(CS)DEFINT(D6,X,0,INFJQ@ 
2 2Pl 
3 SQRTC(3) 


Now let us turn to the problem of finding the poles of 
RCZ). The SOLVE [16] program in MACSYMA knows how to solve 
a number of types of equations. This routine Is used tn 
obtaining the locatfons and multiplicities of the zeros of. 
Q(x), the denominator of R(X). The SOLVE program does this 
by factoring Q(X) over the [Integers and applying formulas to 


each Irreducible factor of degree less than 5. For factors 
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of higher degree only those of the form aX**ntb will be 
solved. This factlilty can be augmented by [Includtng more 
special cases for polynomials of higher degree. If the 
Integrand [s not ratlonal, the problem of findIng poles can 
be much more diffftcult, although SOLVE can find some of 
these also. In any event, the method of residues can be 


successful only If all relevant poles can be located, 


Let p be a root of Q(X) = 0 of multiplicity m which ts 
obtalned employing SOLVE. It Is not necessarily true that p 
Is a pale of R(Z) of order m. This fs because the rational 
function package [10] in MACSYMA does GCD cancellations only 
over the integers while R(Z) ts a rattonal function over the 
complex numbers. Therefore p may not be a pole of R(Z) or p 
may be a pole of R(Z) of order less than m. Fortunately 
WANDERER can pretend that p Ts really a pole of order m and 
proceed with the algortthm for computtIng residues which will 
produce a 0, {ff p ts not a pole, and the correct residue, In 
case p Is a pole of order less than m. More discussfon on 


the computation of residues can be found [tn Sect. 6-1. 


Thus, as the reader can easily vertfy, (D2) of example 
(C2), In page 59, Is valld even Tf A=0 and B=1 (or A=0 and 
B=9). As another example, correct answers were obtained In 


the following integrals even If %1, the square root of -1l, 
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is treated as a simple pole of the fntegrand. 
Ce 

2 
(X = 2 Z1)0X + 1) 
(C7) DEFINT(D6,X,MINF, INF) 
(D7) 00 ee 
(C8) DEFINTC06,X,0, INF) 
(D8) een ennnn 


A spectal case check IS provided as an auxtlfary to the 


above algortthm. 


2 Y dx 
N M 
o° (AX + B) 
P/N - P/N +#M-= 1 SPI 
S(BJAY° AD eee eceemaecee 
M-1 M 
B N SINC3PI P/N) 


for M, N, P positive Integers, M > P/N, P not divistIble by N 
and AB > 0. 
For the case M=1, AB<0, the following Cauchy prtncipal 


value fs used, 
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A 1- P/N &PI ZPl P 
=(- -) (---) COT(----- ) 
B AN N 
This special case Is not necessary for the algorithm 
but ft Is very helpful as far as program efficlency Is 


concerned, especially for large M or N. 


Sometimes substantial computing time can be saved by 
application of differentiatton techniques [fn Integration 


problems. Consider the Infinite Integral 


me dX 
(5) , el ee x 
I 2 10 


o (xX +X + K) 
lf the algorithm for rattonal functions were used straight- 
forwardly, differentiation would have to be carrfed out 18 
tImes to obtain the sum of resIdues tn case both poles lle 


above the real axfs. But since 


jt Is only necessary to differentiate a quadratic expression 


once to obtain the sum of residues to evaluate the Integral 
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ad dx 
2 
ey.) X 


+X + kK 
After thls Its obtalned, 9 more differentiattons with respect 
to a parameter K are needed. I/t means a saving of at least 


8 differentlations. 


In fact, whenever the denominator Q(X) of the !ntegrand 
Is of the form 
v(X)**n 
with v(X) a polynomial and n = 2, this method of diffe- 
rentlatfon may be applicable If the degree of Q(X) Its large 


compared to that of the numerator P(X). 
We have Included this technique tn WANDERER. Let 


n = deg(P(X)), 
s = deg(v(X)), 
QCX) = v(X)*"m, 


Here Is a brief description of the algor!thm. 
'f 2 2mes-n, this algortthm ts not applicable, 
otherwise (m*s-n > 2) 
(1) Set r to the least [Integer = (nt2)/s 
(ti) lfmor 


(a) Compute the ortgiInal [tntegral with m 
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replaced by r and v(X) replaced by v(X) + ZP 
where ZP Is a new parameter Introduced Into the 
problem. Let the answer thus obtained be 
ANS(ZP). 


(b) Return the following as the answer 


m-r 


m-r (r-1)! d 
(-1) 


ma ANS(ZP) 


(m-1)! dZP ZP=0 


(ji) otherwise (m=r), the algorithm Is not necessary 


and thus wlll not be applied, 


The parameter ZP used In thls algorithm makes the 
symbolle differentiation In step (If-b) possIble at all 
times. The final answer Is obtalned by setting ZP to zero. 
Such a parameter wlll be used again later and will be 


referred to as the zero parameter, 


For example, 


1 
(D9) i os a 
2 3 
(X +X + 1) 
(C1LO)DEFINT(D9,X,0,INF)@ 
k ZPI 1 


9 SQRT(3) 2 


In thts example, we have n=0, m=3 and s=2. Thus r=1 


and the value 


Is computed 


ot 
of 


ee 
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2. Infltnite Integrals Involving Algebratc Irratfonal Func- 


tlons 
2.1 Evaluation by Restdues 
Let R(X) be a ratftonal function. Constder 


oo 
D =f X*#*KeR(X) dX, “i K <1 
° 


LIMITCX**(K+1) R(X), X, 0, PLUS) = 0, 
LIMITCX*#(K+1) R(X), X, INF) = 0, 
and R(X) has no poles of order greater than one on the 
positive real axis, then D can be evaluated by applying 


residue theory to the contour Integral 


[ F(Z) dZ, F(Z) = (=Z)¥*K R(Z) 
'e} 


around the contour C shown tn fig. 2. The limit condittons 
given above are convergence tests for D. They are computed 
by DELIMITER. tn computing these Integrals, the Integratton 
program will flrst obtain two sums of restdues S1 and S2 of 


the functfon f(Z) by executing the followtng two steps. 


If R(Z) has poles off the positive real axis, set $l 


to )'’Res(f(Z)) at these poles, otherwise set S1 to 
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0. 
If R(Z) has stmple posttive poles, set $2 to 


yRes(#(Z)) at these poles, otherwise set S2 to 0. 
Now the answer to the Integral D Its given by 


“Pl $1 
steetatetertenteeteateteteteel - ZPl COTC(ZPI! (K + 1)) S2 
SINCZPI CK + 1)) 
For example, MACSYMA produced the following result. 


(C1) 1/CC1+X)eXe%(1/2))@ 


1 
(D1) ee ee ee eee 
SQRT(X) (X + 1) 
(C2) DEFINT(Z,X,0,1NF)@ [3] 
(D2) 3PI 


2.2 Integrals Related to the Beta Function 


From the definttfion of beta functton 
4 
BETACK,S) af X*e(K-1) (1-X)**(S=1) dX, 
) 


where K > 0 and S > 0, One can deduce the following rela- 


tlon 
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where CD > 0 and S$ > K>DO. 


By a further change of vartable tn thls formula, the 
following very useful relatton Is obtatned. 


X BETACA, B) 
fay ... SSaeR Ss dk, Sea Seoss= P) r> 0 


where A = K/r > 0, and B # S-A > 0. 

This rather general formula covers many Interesting 
Infinite Integrals Involving algebrale Itrratlonal functions. 
This ts built Into the programs by usfng special purpose 
pattern recognitton routines (see chapter 6) to examine 
whether the integrand ts of the particular form (a). A 
simplification routine for Beta function was also needed to 


put the results In a s{mpler and better lookIng form. 


(D3) ee 


(C4) DEFINT(D3,X,0,INF)@ [12] 


1S THE EXPRESSION 
K +1 
POSITIVE, NEGATIVE, OR ZERO 
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POSITIVE@ 

IS THE EXPRESSION 

POSITIVE, NEGATIVE, OR ZERO 
POSITIVEG@ 


K 
(Dk) 3 BETACK + 1, = K) 


2.3 Evaluation by Trlgonometric Substitution 
If RCX,Y) Is a ratfonal functtons [fn the two vartables 
X and Y, then an Integral of the form 


a 
U « [nex SQRT(X**2-Axw2)) dX 
3 


can be transformed to an Inftnite Integral of a rational 
function which has been discussed tn Sect. 4-1. The trans- 


formatfon below Is easy to verify. 


Let 
2 2 
SQRT(X - A) 
Y 3 wena nnn nmann 
X+A 
U becomes 
L A(Y +1) 2AY Y dY 
“oA R( ---ne-e-e- eee ee : 
2 2 22 
° 1- Y 1-Y (l=-yY) 
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Now let Z = Y/(1 - Y) In the above Integral, {tt can be shown 


that U Is equivalent to 


Although the above Is all that Is needed for the given 
Integral U, It fs sometimes more efffcftent to transform U 
into a finite [ntegral of a rational function of the trigo- 
nometric functions. That fs, by setting COS(t) = A/X In U, 


one may deduce that 


we A SINC t) 
(iT) U=A R(------ JA TAN(t))  (-"----- ) dt 
) 


cos (t) 
Both of these methods (1) and (If) have been Included 
In the programs. Our heurfstic rule for applying these 


transformations Is as follows. 


Apply (tt) tf it transforms the given Itntegrand [nto 


the form 


COS(t)**#m SINCt)**n, 


otherwise apply (1). 


For example, 
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(D5) —_ mene 
X SQRT(X - A) 


(C6)DEF INT(D5,X,0,1NF)Q 
1 
(D6) = 


2 
A 


2.4 ODifferentlation with Respect to a Parameter 


We have Included two formulas In WANDERER for infinite 
Integrals Involving algebraic trrattonal functlons. They 


are 


7 dX 1 
eee eee ween ne Bo eee ewe meee een wnnnn--= = (A,B,C), 
2 3/2  SQRT(C) (B/2+SQRTCAC)) 
° (AX + BX 4C) 


with A 20, C > 0, and B > “SQRTCAC), and 
e xX dX 1 
Se ee eS SS SS ee R season asens eee tonneeee = V(A,B,C) 
2 3/2 SQRTCA) (B/2#SQRTCAC)) 
° (AX + BX #C) 


where A > 0, C =O and B > = SQRTCAC). 


They are Included because any Integral of the form 
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@ X dX 
(ttt) [ semen anatase aman 
2 N+3/2 
Oo (AX + BX +#C) 


can be computed from these known results. Our tool In dotng 
this ts ditfferentftatton with respect to an appropriate 


parameter. Let Z1, Z2 and Z3 be three zero parameters and 


N SQRTCZPI)/2 
H = (-1) ee, ee ee ee ee . 


3 
GAMMA(= + N) 
2 


The algorithm is as follows. 


Upon dectding that the {tntegrand Is as given tn (lft), 
WANDERER has obtatned the values of M, N, A, B and C. The 
answer to the given Integral Its then computed ustng the 


stmple procedure: 


lf M = 0 and N 0, return U, 


lf M = 1 and N 0, return V, 


If 2N+2 SM, the given Integral Its divergent. 


lf N=M, return 


H ( w-enn--n--=- UCA,B+Z2,C+Z3) 


d 
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If M = Ntl1, return 


N 
d 

H [---- V(A,B+Z2,C), 
dz2 


If 2(N+1) > M > Nt1, then 
(i) Tf M fs even set r = M/2, return the answer 
a aca at a UCA+Z1,B,C+Z3) 
, 
(fl) Tf M ts odd, set r = (M-1)/2 and return 
SSeS a Sais VCA+Z1,B,C+Z3) 


(D7) eee eee 


2 2 
(X +X #1) 


(C8) DEFINT(D7,X,0,1NF)Q 


(D8) wee ne 
6567561 
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3. Infinite Integrals Involving Trigonometric Functlfons. 
3.1 Integrals from Minus Infinity to Infinity 


Let R(X) be a ratftonal functfion which has no real poles 


and 
LIMITCR(X),X, INF) = 0. 


The Integral 


i (41 m X) 
Le ZE R(X) dX, m> 0 
-@ 


Is convergent. Its value can be obtatned by evaluating the 


contour Integral 


(Zio m Z) 
Jj = 2E R(Z) dz 
ra} 


around the contour C given In flg. 1. 


Let Cr be a circular arc with center at Z = 0, radtus r 
and argument t, tl < t < t2. Jordan's Lemma [7] shows that 


If f(Z) approaches 0 unlformly on Cr as r approaches INF, 


(Zl m Z) 
LIMIT ZE F(Z) dZ=0, for tl 2Q, t2% %Pt 
TT —* +00 ¢ 

r 


and that 
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(<1 m Z) 
LIMIT ZE f(Z) dZ=0, for tl 2&Pl, t2 = 2%Pl. 
Y' +00 ¢y 


Therefore L = J as r tends to Infinity and the method 
of residues can be appllted. Moreover, for m < 0 one may 
use the same method by ustng a contour similar to C In the 


lower half complex Z-plane. 


For realm, n and p, let T(X) be SIN(mX), COS(nX), 
9Eee(%l p X) or a functfon Involving sums and/or products of 
these functlons. In complex exponential form T(X) Is a sum 
of constant multiples of functions of the form %E**(Z! k X), 
k real. 


Therefore 


oo 
kK R(X) dX 
hd 


can be Integrated usIng complex contour fntegral. In dotng 
such a problem, terms In the Integrand are sorted into two 
parts. One part requires a contour In the upper-half Z- 


plane, the other a contour In the lower-half. 


Before discusston of additional methods, let us see a 


few examples computed using methods discussed so far. 
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cos(x) 
(D1) 3 
2 
X #1 
(C2) DEFINT(D1,X,MINF,INF)Q (12) 
SPI 
(D2) cas 
%E 
(C3) X*SINCX)/(X**24+1)@ 
X SINCX) 
CDF 2 et ne er 
2 
X +1 
(Ch) DEFINT(D3,X,MINF,INF)@ (12) 
&PI 
(Dk) ae 
SE 
(C5)X*COS(X)/(X**241)@ 
X COS(X) 
COR ean a eI 
2 
X +1 
(C6) DEFINT(D5,X,MINF,INF)@ 
(D6) 0 
(C7) X*COS(X)/(X**24+X+1)G 
X cOS(X) 
COR 2 Fa Reenter eS 
X +Xel1 


(C8) DEFINT(D7,X,MINF,INF)G@ 


(D8) 
1 (SQRT(3) = 31)/2 
(ZP1/2) ((- anne - - %1) ZE 
SQRT(3) 
1 C= SQRT(3)=- %19/2 
+g oS Seesnce ) 2E ) 
SQRT(3) 
(C9) 1/(YE**( SEX) e(X*e2e1) IG 
1 
(D9) et et ee - 
1X 2 
ZE (X + 1) 
(C10) DEFINT (D9,X,MINF,INF)Q 
SPI 
(D10) Soe 
ZE 
(C11) SINCX)*D9@ 
SINCX) 
COLE a i tr 
eI xX 2 
ZE (X + 1) 
(C12) DEFINT(D11,X,MINF, INF)P . 
eP| 
(D12) Beha 
%E 
(C13) D11¥(YEe*(-%1%X)) 
SINCX) 
(D13) eee ween 
2 21 X 2 
%E (X + 1) 


(C14) DEFINT(D13,X,MINF,INF)@ 
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(D14) www 


3.2 Integrals from 0 to Inftntty 


Let us now discuss some tIntegrals tnvolving trigonome- 
tric functions over the range (0 INF), First consider the 


Integrals 


oo n 
tl -f COS(k X ) dX 


9 

aw n 
sek 12 -[ SINCK X ) dX 

co) 


where k a nonzero real constant andn > 1. 


Here agatn, the resIdue theory can be applied to 
evaluate the Integrals. This time the shape of the contour 
Is sllghtly different. The contour Is a sector of a circle 
consisting of a portion of a clrcular arc with center at the 
ortgIn and two stralght lInes jotning the end polnts to the 
origin, (fig. 3) The sector angle depends on n and Is 


“PI/(2n). 


Without loss of generality, let us assume that k > 0. 


Therefore, taking real and Imaginary parts of the above 


relation, we have 
11 = COS(ZPI/2n) G 
and 12 = SINCZPI/2n) G. 


For example 


7/3 
(D15) cos(9 X ) 


(C16) DEFINT(D15,X,0,INF)Q 
3 3 3PI 
3 GAMMA(-) COS(----- ) 
7 
(D16) Re teen enn e nen nn eee 


(C17) SINC 9*X**(7/3))@ 


7/3 
(D17) SINC9 X ) 


(C18) DEFINT(D17,X,0,INF)G@ 


3 3 %Pl 
3 GAMMA(=-) SIN(----- ) 
7 14 


(D18) tccesueaweaeotweaatuee 
3/7 
7*9 


Inctdentally, if ns 1, !1 and !2 do not converge. 
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A more general Integral than {1 and 12 is 


a n 
v=] X EXP(ZI kX ) dX 
° 


where n >--0, Rl(m) > <1, k real and nonzero and n-RIi(m) > 1. 
We shall constder the case k <¢ 0. The case k > 0 Is 


entirely analogous. Let us take a sector-shaped contour fn 


the fourth quadrant as shown In fig. 4. 
y 
Ce x 


cR 


fig 4 


By Cauchy Integral theorem, 


n 
m (Zt k Z) 
Z ZE dZ = 0 
CG 


whtch fmplles that 
fos) 


m n 
J = EXP(-4P I aneaess/cann| R EXP(k R ) dR 
° 


Ss 
= EXP($3PL 31 (m+1)/(2n)) GAMMA(------- ) 


s 
n(=-k) 
where s = (mtl)/n. 
Let us look at an example. 
(2 31 + 3) 
xX 
(D19) ee ee eee 
3 
(31 X / 2) 
ZE 
(C20) DEFINT(D19,X,0,INF)@ 
2 3t + & 
GAMMA(-------- ) 
3 
COZO)) © earner eas seca eevee ria ea eae 
(2 231 +h) / 3 
3 (%1 7 2) 
A stmllar result for k > 0 ts given by 
(m+1)/n 
J= EXP(ZPL ZI (mt1)/(2n)) GAMMA(----------- ) 
(m+1)/n 
k n 


The following [Integrals can be obtatned readily from 


the above relatton. 


aed 00 
n m n m 
[cose xX )X dX and [sine x )X dx 
° 
0 


For Instance, 
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aaa g “Ss r(r- 1) (? re -(s-2) 
SIN (X) X UX = wceweccene SIN (X) X dX 
_ (s-1)(s-2) Jo 


r @ r ~(s-2) 
¢ scam tne: Eom SIN (X) x dx, 
(s-1l)(s-2)Jo 


where r > (s"1) > 1. As one can see, repeated application 
of this relation wlll reduce U to a sum of Integrals of the 
same form as U but with N = 1 or 2. For N#1, AUDITOR uses 


the following formula 


co U/Z 
p “1 p-l 
[sm (Xx) xX dX -[ SIN (X) dX, p > 0 odd. 
oOo 


° 


The tntegral on the right hand side can be evaluated eastly. 
Integrals of this type are considered fn Sect. 5-2.2. For 


N=2, we have 


q 2% ZPI1 fq-3/2 
SIN (X) X dX = --- 
oO 2 q-1 


where q ts an Integer 22. For example 


SINCR X) 
(DI) en 


88 


(C2) DEFINT(%,2%,6, INF)@ 


2 | 
(D2) ae 
2 
(C3) SINC T+X)#a2/(Xe¥2)Q 
2 
SIN (Q X) 
(03) enn 
2 
x 
(Ch) DEFINT(%,X,0,INF)G@ 
Pl Q 


OS 
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4h. Inftnite tntegral Involving Logartthm Functltons 
4.1 Evaluation by Recursfon and Contour Integratton 


Let R(X) be a rational functton whIch Is even fn X. 


The Integral 
2 
VON) -f LOG (X) R(X) dX, R(X) = PCX)/QCX), 
o 


with deg(Q(X))-deg(P(X)) 2 2, can be evaluated by applying 
residue theory In a recursive manner. Recall that LOG fs an 
abbreviatlon of PLOG when the argument ts real and positive. 


Consftder the contour Integral 
N 
JCN) = F(Z) dZ, F(Z) = PLOG (Z) R(Z) 
d 


where C ts the Indented contour tn flg. 5. As ¥ approaches 
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1) Declare each of AV and AJ to be an 1 XN array 
2) Set |! to 0 then go to step 4 
3) Set AJ(Tt) to JCI)/2 


h) If I N, return the answer V(N) 


5) Set AV(I) to V(t) 


6) Set { to [+l then go to step 3 
Computation of V(1): 


If t = 0, return the value of the Integral 


@ 
[rw dx, 
° 


(storing away [nformatton about poles and 
corresponding residues of R(Z) for computation of J), 


otherwise, compute from the formula 
i 


1 I K K 
V(t) = AJ(T) - - | Pl %t AVCT = K), 
2 K 
K> 
Computatton of J(1)/2 : 


Ustng existing Itnformatton of poles and restdues of 


R(Z) compute by reS8Iidue theory from the formula 


ZPI %1C }° Res(F) over poles of R Inside C), 


For example, 


o2 


2 
LOG.¢X) 
(D1) = seeeeee 
2 
X +1] 
(C2) DEFINT(D1,X,0,INF)@ {1] 
3 
%P | 
(D2) elas 
8 


4.2 Method of Differentiating and Introducing Parameters 


A very usetul method for Integrands [tnvolving LOG(X) fs 


differentiation. Consftder an Integral of the form 


ror) 
K 
ACK) -f R(X) X LOG(X) dX, K # 0 and =1.¢ kK < J, 
° 


where R(X) [Is rattlonal tn X and k ts a parameter which 


occurs nowhere else in the [ntegrand. If the Integral 


@ K 
B(K) i R(X) X dX 
Oo 


Is convergent and can be evaluated, then ACK) ts given by 
d 
ACK) = -= BCK), 
dk 
This method of differentiatlon is valid because B(K) is 
convergent and A(K) Is unlformly convergent for every closed 


interval contalned In the set of polnts (K | K#0 and -1<K<1), 
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For example 


K 
X LOG(X) 
(D3) -------+- 
X +3 
(C4) DEFINT(D3,X,0,1INF)Q [12] 


K K OD 
(D4) LOG(3) 3 BETAC(K+1,- K) + 3) (=-BETACr, - K) 
Dr 
D 
- =--BETACK + 1,r) ) 
Dr r=-K 


r=(1+K) 


In the foregotng dIiscusston, K has been assumed 
symbolic. Yet, tn solving actual problems K may very well 
be a number. This difficulty can be overcome by Introducing 
a zero parameter ZP, We first replace X**k by X*#*(k+#ZP) fn 
F(X). This permits us to proceed as above and then diffe- 
rentlate with respect to ZP, After differentlation the 


result Is then evaluated at ZP = 0. 


The property of the logarithm functton 
~LOG(X) = LOG(1/X) 


can sometimes be used [In evaluating Integrals of the form 


oS 


b 


J -[ LOGCg(X)) FCX) dX 
a 


Is convergent, then 


b > ucxre' (xX) dx 
FUORI Sh ee once sees ), 
adja 2(X) 


The SIN tntegratton program [20] tn MACSYMA Is used to 


compute the Indefinite Integral U(X). WANDERER uses thls 


method to evaluate Integrals of the form 


cient E a 
X LOG(1 + X ) dX, atl > LO 1. 
3° 


Here Is an Integral evaluated by the method of {ntegra- 


tion by parts. 


(D7) eee ee 


(D8) eee n 
5 
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5. Infinite Integrals Involving Exponential Functions 
5.1 Method of SubstIitutton 


Integrals with Integrands which are ratltonal functlons 
of ZE**(K*X), K real and non-zero, are relatively easy to 
compute. Without loss of generality, let us assume K > 0. 


and constider 


@ 
(1) [ eer dx. 


00 


If one makes a change of varlable 


(2) Y = ZE*e(KeX), 


@ 
cur f RCY)/Y dY. 
fe) 


This tntegral converges !f (1) does. That {fs to say 


then (1) becomes 


R(Y) has Y as a factor If (1) converges. For the Integral 


ee) 
(3) [reece dX 
) 


one can make a stmllar change of varlable 
(4) S + 1 = ZExn(KeX) 


which gives 
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oo 
(5) [esenisen ds. 
° 


Integrals of the form of (3) and (5) can be evaluated 
by contour Integral and residue theory as Indicated In sect. 


1. The methods (2) and (4) can In general be applied to any 


b 
[ tec dX 
a 


and the resulting ftntegral may be much simpler. For 


Integral of the form 


Instance 
1 
(D1) a a wat to re ee 
X/3 5 
“Ee (aeec8 + 7) 
X/3 
SE 
(C2) DEFINT(D1,X,0,INF)IQ 
7 
3 LOG(--) 
12 
(D2) iss ees ecw eee 
5 
X/4 
zE 
(D3) eee 
X/2 
9 %E + & 


(Ch) DEFINT(D3,X,MINF,INF)Q 
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(D4) --- 


5.2 Use of Contour Itntegratton 


A somewhat more [Interesting Integral Its 


00 
| f RCZEw*X) P(X) dX 
-09 


where P(X) is a polynomial and R(X) a ratlfonal function with 
complex coefficients such that 
LIMITCROZE**X),X, INF) = 0 
and 
LIMITCROZE**X),X,MINF) = 0. 
Let us flrst determine a polynomial Q(Z) with complex 
coefficients such that 
(6) Q(X) =~ Q(X + 2 SPI %1) = P(X). 
Q(Z) exists and can be computed by the method of 


undetermined coeffictents. Now constder a contour Integral 


J a R(ZE**Z) Q(Z) dZ 
/C 


taken around a rectangular contour as shown [In flg. 6. 


As the absolute value of X approaches [nftntty, 
contributlons from the vertical segments of the rectangular 


contour vanish. Hence, 


co co 
J -[ RCZE**X) QCX) dX -[ RCZE**X) QUX+2 2PL 41) dX, 


CO 00 


It follows from (6) that l=J. Now the problem of evaluating 
the integral | has been reduced to finding the residues of 
RCZE*x*Z)Q(Z) for 0 <= Im(Z) < 2 YPI. To do thts, the poles 
of R(Z) are obtatned flrst. If w ts such a pole then 
GLOG(w) 

where GLOG stands for the branch of LOG with Imaginary part 
between 0 and 2*%P!, Is a pole of the same order for 
RCZE**Z)*Q(Z) InsIde the closed contour, Ajll such poles can 


be obtalned tn this manner. One example ts 


(D5) ene 
SINH(X) = 31 


(C6) DEFINT(D5,X,MINF,INF)Q 
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(D6) %PI 


Note that the Integrand (D5) has a pole of order two at 


Z= ZPIl ZI/2. 
5.3 The GAMMA Functton and Related Integrals 


A very Important functlon closely related to the 
evaluation of Infinite [ntegrals [Involving exponentials Its 


the GAMMA functlon generally defined as 


oO 
GAMMA (Z) -[ ZEwe(-t) tee(Z-1) dt, R1(Z) > 0. 
o 


Also of use {ts tts logarithmic derivative, the PSI 


function 


d 
PSI(Z) =(€=- GAMMA(Z))/GAMMA(Z). 
dz 


A stmplltficatton routtIne for GAMMA functton has been 


written to make use of the many properttes of this functton. 


From the definition of the GAMMA functlon, one can 


derive the followIng very useful relatton 
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where a = (D + 1)/B, RICA) > 0, RICD) > -1 and R1(B) 
nonzero. 

WANDERER has programs deslgned to recognize this form 
and return the result. Of course these programs have to 
examine the stgns of the relevant quantities carefully 


before generating an answer. 


For example: 


-X 
(D7) 43 


(C8) DEFINT(D7,X,0,1NF)G@ 
SQRT(ZP1) 


(D8) eonennn- 
2 


5.4 Integral Related to the Laplace Transform 


Let f(t) be a function of a real vartable t, then [ts 


Laplace transform L(f(t)) Its defined as 


@ 
(=t s) 
LCf#(t)) = F(s) =| SE f(t) dt, 


° 
Many of such Integrals can be evaluated by the pro- 


grams. Here are some examples produced by WANDERER. 
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CHAPTER V 
INTEGRALS OVER A FINITE RANGE 


0. Introduction 


Many finite Integrals are proper Integrals whose 
Indefinite Integrals exist In closed form and can be com- 
puted rather easily. For such an ftntegral, the evaluation 
method Is very stralght-forward. WANDERER stImply obtalns 
the corresponding Indefinite Integrals and then substitutes 
the limits of Integration. The antiderlvatives are computed 
by use of SIN. This method will be referred to as the antI- 


derfvattve method. 


A finite Integral 


b 
[ dX 
a 


Is Improper If f(X) becomes Infinfte at some poltnt c, 

b 2c 2a. In order to avold having to spend computation 
time looking for an unknown number of sIngulartties of f 
between a and b, WANDERER wlll assume that Integrands of 


finite Integrals, with the exception of rational Integrands, 
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can become Infinite only at the end points of the range of 
Integration. This ts not a severe restriction and does not 
decrease the number of Integrals [t can handle. This ts 
true In the sense that any glven range can be subdivided to 
conform to the above conventlton. This conventfon makes [t 
easy to determine whether a given finite Integral Is 
improper. WANDERER stmply checks the value of the Integrand 
at the limits of Integration. If the given Integral Is 
Improper, Its divergence Is tested before any attempt at 
evaluation. WANDERER uses a limit test for absolute 
divergence which Is discussed In Scet. 6-6. If the given 
Integral diverges, WANDERER will so Indicate In the output. 
If the antidertvatIve can be computed, then the answer Is 
sometimes obtalned by employing a limiting process when 


substituting the upper and lower limits of Integrattfon, 


In thts chapter, attentton will be focused on the 
definite Integrals whose corresponding Indefinite Integrals 


are difficult to compute or do not exist. 


1. Finite Integrals of Rational Functtons. 


For Integrals such as 


b 
U -[ R(X) dX 
a 


where R(X) Is rattonal In X, a and b are finite, WANDERER 
computes U by transforming [ft Into an Infinite Integral by a 


change of vartable. 


Let us write 


(X = a) (a + Y) 
(1) Y = wenwnee ° X = ween nwee ? 
(b = X) (1 + Y) 
(b - a) dY 
dX & cceeenence 
(1 + Y) 


2 
(Y + 1) 
This Integral can be tntegrated readily by methods of con- 


tour Integration and other means discussed In sect. 4-1. 
As an example let R(X) be the expresston (D1) 


2 
x = 3 
The Indefinite Integral of R(X) from 0 to 1 computed by 
the substitution method above fs given by 


(C2) DEFINT(D1,X,0,1)@ 


LOG(2 = SQRT(3)) 
(D2) et tt nee nen 
2 SQRT(3) 
One knows that the Indefitntte Integral of (D1) extsts. 
In fact the followtng has been obtalned ustng the command 
INTEGRATE fn MACSYMA 


(C3) INTEGRATE(D1,X)@ 


(D3) (LOG(2 X = 2 SQRT(3)) - LOG(2 X + 2 SQRT(3))) 
/(2 SQRT(3)) 
The reader may easlly obtafn a result equivalent to 


(D2) by Substituting In (D3) the limits of Itntegratfon,. 


2. Ratlonal Functtons of Trigonometric Functions 
2.1 <A Typical Applicatton of Contour Integration 


If RCX,Y) ts a rattonal function [tn two varfables X and 


Y, an tntegral tn the form 
2K 
U = [ Rccoscxy,sincxr) dX 
°o 


is eastly transformed to an [tntegral around a closed con- 
tour. By setting 


2 2 
zZ +1 Zz -1 
COS(X) = ------ ; SIN(X) = ------ , 
2Z 2412 
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dZ 
and dX = --=-, 
%1 Z 

the glven Integral U becomes 


2 2 
1 Z+1 Z-1 dz 
wn | F(wneene poaeeee )} =. 
tiaj 22 23tZz 2 


where Cl is the pos!ttIively orfentated unfit circle with 
center at Z = 0. This contour [tntegral can then be 


evaluated by finding the sum of res!Idues Instftde the circle. 


Actually this transformatton can be applited [n general 


to any Integral In the form 
2n 
f rceeeecet X)) dX 
oe 
where R Is rational, by the change of varlable 
Z = ZEwwe(Zl xX) 


The transformatton process Its simple. The key pofnt tn 
this algorithm Is the determinatton of whether the Integrand 


is In fact a rational functlon of ZExx(3!l xX), 


In calculating the sum of the restdues, only poles 
Instde the contour contribute. SImple poles on the unit 
cIrcle, t.e. those with absolute value 1, cause the 


princtpal value to be computed. The Integral fs divergent 
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If poles of order greater than 1 happen to be on the circle. 


Examples: 
(Ch) COS(X)*#2-SIN(X)Q 
(D4) cos (x) - SINCX) 
(C5) DEFINT(DH,X,0,2*%P1)@ 
(D5) ZP1 


(C6) SER 2eZi eX) /(ZEwe (ZL *X)+3)Q 
2 41 X 

(06) eee nnn 

(C7) DEFINT(D6,X,0,2*ZP1)@ 


(07) 2 2PI 


2.2 Uttlizatiton of the PertodictIty of the Tri gonometric 


Functions 


In this sectfton Integrals of functfons Involving trigono- 
metric functlons over a variety of ranges will be 


considered, 


Let T be a functton of X deftned by 


M N 
TCX) = COS (X) SIN (X) 
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where R1(M) > -1 and RI(N) > <1. The following formula can 


be deduced from the defitnitton of Beta function. (Sect. 4- 


M + + 
js/2 GAMMA (¢ mee ) GAMMA( —e ) 
(1) TOX) dX = sanccn neces e ce cw ewes eneea 
3 2 GAMMAC(N + M + 2)/2) 
A stmplte example Is 
1/3 1/2 
(D8) cos(Xx) SIN (X) 


(C9) DEFINT(D8,X,0,%P1/2)@ 


2 3 
6 GAMMA(-) GAMMA( =) 
3 4 
(D9) a art ar at am ° 
5 
5 GAMMA(-=-) 
12 


The usefulness of (1) Is Increased by the fact that It 
Is posstble to express definite Integrals of T(X) over a 
vartety of ranges tn terms of that of T(X) over (0 3PI/2). 
For Instance, the following relations are true for any func- 


tion f. 
t 
il F(SINCX),COS(X)) dX = 
Oo 


WZ 
[ crcsincas,coscan + FCSINCX),-COS(X))) dX; 


(o] 


3/2 
[ tesinc,c0s009 dX = 


R R/2 
J ecsimcas e080) a [ f(-SINCX),-COS(X)) 
ro) ° 


More generally, let f(X) be a pertodic function of X 
with perlod 2%P!. That Its 
F(X + 22PI1) = F(X), 


An integral of f(X) over some range (a b), 


b 
S 5 ‘| F(X) dX, 
a 


can always be written as a sum of ftntegrals In the form 


2 


c d 
(2) n f F(X) dX “f f(X) dX oi) £(X) dX 
fo] o 


° 


for some Integer n and 2%P1 > c, d2O. This ts true for 


there exist [Integers p and q such that 


a= 2p %Pi +d, 


and b 


2q PI + c. 


Then S Is equivalent to the sum (2) with n = (q - p). 


Programs have been written to perform this reduction 
and they are applied when the Integrand has a pertod 2%Pl 


and the difference (a - b) has ZPI as a factor. some 
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dX, 
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examples computed by WANDERER are 
1/3 2 
(D10) cOS(X) SIN (X) 


(C11) DEFINT(D10,X,-%P1/2,%P1/2)@ 


2 
18 SQRT(ZPI) GAMMA(-) 
3 
(D11) ee eee een 
1 
7 GAMMA(=) 
6 


(C12) COS(X)**35*SIN(X) #*#2Q 


3 2 
(D12) COS (X) SIN (X) 
(C13) DEFINTC(D12,X,3*%P1/2,3*ZPI1)@ 
2 
(D13) == 


15 
3. Finite Integrals of Algebraltc !rrattonal Functtons 
3.1 Ratlonaltzting the Integrand 


If RCX,Y) ts a rattonal function In X and Y, the Inte- 


eral 


. 2 2 
(1) K -[ R(SQRT(A = X ),X) dX, 
b 


where A =>c > b 3-A, can be ratlonallzed. 


Let us write 


2 2 
A = SQRT(A - X ) 
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(2) Y 3 scene enn ne en nee = v(X), 


which gives 


(3) Xs secee- 


SubstitutIng (3) for X, (1) becomes 


cb> 


which Is an Integral of a ratfonal function. 


thls type have been discussed In Sect. 1. 


Stmilarly for the Integral 


¢ 

2 2 

J | R(X,SQRT(X = A )) dX 
b 


where c > b 2A > OQ, the change of vartable 


can be made to convert J to 


2 2 
[" ACL -Y) 2AY (1-¥ ) dY 
wi 


2 2 2 2 
Ee ame 1+ yY (1 +yY) 


Integrals of 
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which ts also an Integral of a rational functlton. Here are 


two such Integrals evaluated by WANDERER, 


1 
(D1) teen nnn 
2 

X SQRT(X = 9) 

(C2) DEFINT(D1,X,3,4)@ [12] 
3! 3 + 31 SQRT(7) 
(D2) e-- PLOG(------------- r) 
3 4 


(C3) 1/CCX+1)*Ch-X¥e2)%%(1/2))8 


(D3) ee ee ene en nn 
2 
(X + 1) SQRT(H - X ) 


(Ch) DEFINT(D3,X,0,2)@ 


LOG(2 + SQRT(3)) 
(Dh) a a a ee ee a ae ee ee 
SQRT(3) 
Another method Is to try to transform the given Inte- 


gral to an Inftnite Integral by the change of vartable given 


In (1) of Sect. 1. For Instance the tntegral 


B 
[cx sonrces ~ A)(B - X))) dx 
A 


can be converted to 


The method of rattonalfzatton can also be applied to 


b 
[ro (CX + D)**(1/Q)) dX 
a 


where Q [Is an Integer, C and D are constants.and the range 


of integration needs not be finite. 


For Integrals of this type the substitution 


1/Q 
Y = (C X + D) , 


will convert the given [Integral Into that of a ratfonal 


function tn Y. 
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The princtpal task In thls converston method fs to 
recognize that the given Integrand Its of the part!icular form 
and to obtain the Integer Q. In essence, the algorithm for 


dotng this fs 


(ft) Obtaln a llst L of all distinct trrattonal 
parts In the glven Integrand, 

(it) If elements In L are fractional powers of an 
Identical lftnear polynomlal In X, the pattern ts 
matched and Q Is set to the lcm of the denomInators 
of all the exponents of the linear polynomtal, 


otherwlse the pattern [s not matched. 


3.2 Integrals Related to the BETA Functlon 


The BETA functlon Is deflned by the Integral 


4 ket L+2 
BETACK,L) = Xx (1 - X ) dx 


0 
for RICK) > 0 and RICR) > 0. From this definition one may 


readily deduce the relation 


. 1, - 1 cLe-1 1 K 
(3) X (1 - xX) dX = = BETAC-,L), 
3 Cc Cc 


WANDERER appltes this formula by recognizing the form 


of the given Integrand. [nn chapter 6 some techniques and 
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programs for pattern recognition are discussed. For inte- 


grals In the form 


B 
aoe | C Coke 
(X = A) (B =X) dx, 
A 


a simply substlIitutton, 


Y = (X - A)/B, 


wlll transform it [nto (3). 


4h, Fintte Integrals tnvolving Logarithm Functtons 


When the Integrand ts a functton of LOG(X), a given 
Integral may, In many cases, be evaluated by transforming [t 


into one which Involves exponenttal functftons. ConsIder 


b 
U gi f(LOG(X)) dX, b>a 2Q, 
a 
the substitution 
~Y 
(1) X = ZE for 1 >a, 
or 
Y 
(Tt) X = ZE fora 21, 
converts U to 
s -Y 
fren ZE =e, r= -LOG(b), s= -LOG(a), 
tT 
or 
s Y 
fron ZE dY, r=LOG(a), s= LOG(b), 
T 


respectively. In case a = 0 and b = 1 (or INF), the use of 
C(t) or CPt) will result tn an tnfltnite Integral which can 
often be evaluated readily using methods provided in 


WANDERER. For Instance by use of (1) the [Integral 


$  q 
f rocex ) dX 
o 
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becomes 


(o ¢] 
q -Y 
[ow ZE dY, 
° 


which ts an tnfinitte Integral already studied (Sect. 4-5.3), 


The method (f) or (tf) can also be applied to tntegrals of 


b 
‘ 
[: F(LOG(X)) dX, 
a 


4 re) 
Sl R S$ -RY-e=-1 
LOG (=) X dX = Y $&€E dy, 
Po) x ° 


For Integrands ftnvolving LOG(F(X)), It may sometimes be 


the form 


as In 


possible to simplify the given Integral by the substitution 
Y= €(X), X = gtY) 
where the function g Is the Inverse function of f. A very 


simple application of this method Is conversion of 


[e+] 
; x*LOG(SQRT(X)+a) dX 
a 


2 
ad 2 
af (Y-a) LOGCY) dY. 
oO 


Techniques of differentiation with respect to a 


to the integral 


parameter, as detailed th Sect. 4-4.2, can be employed for 


finite Integrals with a factor 


K 
X LOG(X), 


In the integrand. 


(D1) 
(C2) 


(D2) 


(C3) 


(D3) 


(Ch) 


(Dh) 


(C5) 


(D5) 


(C6) 


(D6) 


Examples 3: 


K 
LOG (X) 
DEFINT(D1,X,0,1)@ 


K 
(-1) GAMMACK + 1) 


LOG(X) **(1/2)/(X**2)@ 
SQRT(LOG(X)) 


DEFINT(D3,X,1,1INF)G@ 
SQRT(ZP1) 


1/3 
X 


ee ee ee ee ee ee er se 


SQRT(-1) SQRTCLOG(X)) 


DEFINT(D5,X,0,1)@ 
SQRTCZPI) 


RICK) # 0, 
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[12] 


[12] 
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(C7) LOGCX}*t Cince etl s2))/X) ** (1/2) 


SUPTCL = -SQRECK)) LOGCX) 


(O70) - q - #ergeehiemeseercsers se 
SQRTCX) 
(C8) DEFINT{O/7,«£,0,1)4 
D 3 i 
(03) u (=-BETACr,-)) 


Dr 2 rel 
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(1) Set V to the quottent of V(Z)/(Z = p)**m 
whtch Is computed by long diviston. 


(if) Return the result computed from 


Z=p 
otherwise, (m > 1, V not a polynomial) compute and 
return as answer the limlt 


m-1 


1 d m 
Lin ———— (| (Z = p) f(Z) 
Zp (m- 1)! Id 


1.2 Evaluation of Contour Integrals by Residue Theory 


In order to evaluate an Integral of f(Z) around a 
closed contour C by residue theory, {ft Is necessary to 
locate all poles of f Inside C. After this Its done, the 
remaining problem Is to compute the sum of residues of f at 
these polInts efficiently. The difflcult part Is finding 
poles. WANDERER employs the SOLVE routine In MACSYMA_ to 
solve V(Z)=0. 


For V(Z) a polynomlal In Z, SOLVE finds Its zeros by 
factoring over the Integers and applying formulas to each 
factor of degree less than 5. For factors of higher degree 


only those of the form a*Z**ntb will be solved. The problem 
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of factoring large arbitrary polynomlals [Is non-trivial, to 
say the least. The development of a more powerful factoring 
algorithm which will factor polynomials over a larger ring 
than the Integers would certainly be helpful to SOLVE and 
WANDERER, The location of zeros Its usually more difficult ff 
V(Z) Is not a polynomial In Z but some more comp) tcated 
function. In such a case, WANDERER usually uses methods 
other than the residue theory. An exception Is when V(Z) is 
a polynomial In %&**Z which has been discussed tn Sect. 4- 


5.2 


Thus within the limitations of SOLVE, V(Z) will be 
solved and Its zeros sorted Into a list of pafrs. Each patr 
contalIntng a zero and {ts mult!Ipl{ictity, such as 


L = ( (Zl,m1) , (Z2,m2) , (Z3,m3) , 1.02 


L ts then sorted Into two lists Ll and L2, dlscardtng 


poles outside the closed contour C, such that 


L1 


a list of all stmple poles 


L2 


a lfst of other poles palred with thelr 


multiplictttes. 


At this pofnt, we can apply the RESIDUE ALGORITHM to 
obtain the desftred sum. To avotd repeated calculation of 


V'(Z), the program checks whether L1 Is empty. If L1 Is not 
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an empty lIst, V'(Z) will be computed and stored for 


possibly repeated reference later In the computation. 


Note that It has been assumed from the begInning of 
this section that a pole of f(Z) would be a zero of V(Z). 
This assumption {ts quite reasonable for almost all of the 
applications In evaluation of definite Integrals of 
elementary functlons. However a zero of V(Z) need not be a 


pole of f(Z). For tnstance Z = %! ts a zero of 


but not a pole of 


(Z + $1) 


For each zero p of V(Z) we may check the value of U(p) 
to see If p Is really a pole of f(z). Although It Is not 
clear what can be done If U(p) = 0, since p may still be a 
pole of lower order. A better method Is to Ignore the fact 
that If V(p) = 0, U(p) may also be 0 and pretend that p Is 
an actual pole of f(Z). This ts valld because a residue at 
any removable stngular potnt will turn out to be 0. 


Furthermore, the residue obtained at a pole of order m Is 
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2. Obtaining Real and Imagtnary Parts 


In the course of evaluating a definite Integral by com- 
plex contour Integratlon, the need to take the real or the 
imaginary part of an expression often arise. For Instance, 
to see whether a pole, p, Iles above or below the real axis 
the slgn of Im(p), the Imaginary part of p, ts examined. To 
determtIne If p Illes Insftde the unlit circle at Z=0, It Is 
needed to compute ABS(p) which Involves taking the rea) and 


Imaginary parts of p. 


The algorithm for obtalning R1l(p) Is presented as a 


representative of simllar procedures used. 
Algortthm REALPART(p) : 


1) {tf p ts $I return 0, If p Ts a number or any 
other atomic symbol, return p. 
2) If p ts a sum, (p “> pj) then return 
}. REALPART( pj) 
3) If Si a product (p = pl¥p2) then return 
REALPART( pl) *REALPART( p2) 
- IMPART(Cp1)*IMPART( p2) 
4h) If p = %E**pl return 


SE*x*REALPART(p1)*COSCIMPART(p1)) 
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5) If p = pl**p2, convert p to %E**(p2*LOG(pl)) then 
go to step 4, 
6) If p = LOG(pl), return LOGCABS(p1)) 


7) Otherwise, return the form RI(p). 


3. A Heurlstic Pattern Recognition Program 


It Is often the case that some pattern recognition fs 
needed, at one stage or another durtng the evaluation of an 
Integral. Although many Integrals can be evaluated without 
any pattern recognition, this capability remains Important 
to WANDERER. One specific pattern shall be discussed as a 
representative of such methods tn WANDERER. Constder the 
pattern 


N M 
P(X) = (B X + A) 


where B, N, A and M are free of X (f.e. do not {nvolve X) 
and all except A must be non-zero. This pattern should 
match, for example, every one of the expresstons 
1/2 -1/3 
Xx, 2X +1, (X + 21) 
9) q Z 
(X - 1), and xX #2 X + Ll, 


The last expression Is the expanston of 
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2 
Cio WD 
Expresstons In expanded form present some trouble for 
the recognittion algorithm. Thts dlfflculty Is overcome by 
the use of differentlatlon and rational simplification. 
That ts If an expresston E(X) ts equivalent to an expression 


matching the pattern P(X), we can compute 


E'¢X) / ECX) = UCX) / VOX), 
cancell!Ing all common factors [In the numerator U(X) and 
denominator V(X). Then VCX) should match p(X) with M2 1, 
As a result of thls match, some values are assigned to the 
varlables A and B. By use of these values the correct value 
of M and N can be recovered from U(X). The values for A and 
B thus obtalned ray differ from the true values by a 
constant factor, This happens whenever these true values 
have common factors. The real values of A and B can be 
determined by comparing A**M to E(0). This procedure was 


suggested by Moses. 


The full algorithm used for matching P(X) will be 
described. Let It be called PM, PM uses another routIne PN 


which recognizes the pattern 


N 
P(X) = (B X + A) 


As one can see In the following algortfthm, the values 


at 


Ll/m l/m 
AsAr and B = 8 r, 


then P(X) ts matched. 


Algorithm PNCECX),X) 

1) tf E ts free of X, the pattern not matched, 
2) Match A to £(0) and set E to E- A, 

3) lf E ts of the form reX*xs, the pattern Is 
matched (B to r, N tos). Otherwlse there Is no 


match. 
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4h, <A Procedure for Change of Varfables 


Substitution of a new varlable Y for a subexpression, 
say g(X), of the Integrand In a given Integral ts a 


frequently used method In Integratfon. Let the given Inte- 


b 
J -[ F(X) dX, 
a 


then the transformed Integral would be [In the form 


d 
J -[ F(Y) dy 
c 


where F(Y) might be considerably sfmpler than f(X). FCY), ¢ 


gral be 


and d are computed by a procedure which Is called whenever a 
change of varlable Is needed. tt makes use of two other 


modules of MACSYMA, namely SOLVE and DELIMITER. 


1) Use SOLVE to solve for X tn Y = g(X), obtaltning 
X = h(Y), the Inverse function of g. 

It) tf h can not he obtalned, return Indication of 
fallure, 


titi) Compute, using the DELIMITER, c and d as tn 
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LIMITCZ(X), X, a, PLUS), 


a 
u" 


d = LIMITCg(X), X, b, MINUS), 


Iv) Obtaftn fFCY) by assigning It a value computed 
from 


FChCY)) AMY), 


5. Solving Systems of Linear Algebratic Equattons 


In sect. 4-5.2 the need to compute a polynomial Q(X) 


from a given one P(X) satisfying a given relation 


(1) QCX)-Q(X + 2 SPI 41) = PCX) 


has been mentfoned. 


The method of undetermined coefficlents ts used to 
determine Q(X). Let Q(X) be a polynomfal In X with degree 
one higher than P(X) and unknown coefflclents CO, Cl, ... , 
Cn. That Is 

n+l n 
Q(X) = X + Cn X +... + C1 X + CO. 
Equation (1) will give n lfnear relations among these 
coeffictents. Therefore the value of these C's can be 


obtalned by solving the system of linear equations they must 
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satisfy. 


For the purpose of solving systems of linear algebraltc 
equations, A method known as the "Two-step fractlon-free 
Gausstan elimination" [14] has been Implemented. ThIs 
method Is an Improvement over a correspondtng one-step 
method and features a procedure that keeps the size of the 
Intermedfate expressions In the course of the reductlon down 
by dividing them by a common factor which the procedure can 
predict. The advantage of this method over a more 


effictent and elaborate scheme [18] fs Its sImplicIty. 
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6. Convergence of Integrals 


There are many methods to determine the convergence or 
divergence of a given Improper Integral. Some are limit 
tests, others comparison tests etc. Some test for absolute 
convergence, others conditfonal or uniform convergence. It 
ts a posstble area for future work. There Is no real need 
to have such elaborate schemes tn WANDERER, for after all [tt 


Is the value of the given Integral that ts destred. 


b 
J -[ F(X) dX, 
a 


If f(X) ts a rational functlton of X, WANDERER combInes 


Constder 


convergence tests with evaluation algorithms as explalned [n 
chapters & and 5. If f Is not rattonal, then [t ts not 
allowed to become Infinite except at a and b. This Is a 
convention on Inputs used by WANDERER. Thus, J Is a proper 
integral !f a, b, f(a) and f(b) are finite. If J Is 


Improper, WANDERER has a test for absolute divergence. 


Test for absolute divergence: 


(1) If a and b are finlte and b> a, then J Its 


absolutely divergent If 


LIMIT(C(Cb=X) F(X),X,b,MINUS) # 0, 
or 
LIMITC(X-a) f(X),X,a,PLUS) # 0. 
(2) If ats Inftntte, then J is absolutely divergent tf 


LIMITCX¥f(X),X,a) # 0. 


(3) If b ts Infinite, then J Is absolutely divergent if 
LIMITCX*f(X),X,b) # 0. 
If a functton G(X) exists such that G'(X) = f(X), then 
the method of antiderlvative can be used. That Is to com- 


pute J by evaluating 
LEIMIT(GCX),X,b,MINUS) = LIMITCGOX),X,a,PLUS). 


If this value Is finite, It Is the value of the given 


Improper Integral. If tt does not exist, then J diverges. 
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CHAPTER VII 
AN APPLICATION OF MACSYMA 
AND WANDERER 


0. Introductfon 


The usefulness of a general purpose algebraic manfpula- 
tlon system such as MACSYMA In factlitatIng the solution of 
mathematical problems has been demonstrated. Using his 
"Symbollc Mathematical Laboratory" [15], Martin demonstrated 
solutions to three demanding problems [In applied 
mathematics. These examples emphasize the fact that routtIne 
algebraic computatlon can be done by computer programs not 
only without error but much faster than by hand. The value 
of such a system Is especltally appreciated when the expres-~ 
slons Involved are large and complicated. By employing 
such a computer facility the human problem solver may be 
freed from the tedious and untnsptIring manltpulations to 
think more about the profound aspects of his problem. But 
this ts not all such a system can do. Moses's SIN [20], a 
program for Indefinite Integratton, provides a good example 


of successful mechanizatfton of a mathematical process whlich 
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Is far from routtne or straight-forward. The [ntroductton 
of SIN broadened the scope of algebralic manltpulation systems 


stenificantly. 


The purpose of this chapter ts to show how MACSYMA and 
WANDERER can be used to help solve complftcated problems that 
are of practical tmportance. One such problem Is the asymp- 
totic evaluation of certain contour [Integrals artstIng In 
mathematical physics. Usually one starts with one or a set 
of differential equatftons describing a phystcal problem. 
Solving these equattons by one method or another, most often 
by Integral transforms, one will arrive at a solutfon fn the 
form of a definite Integral whitch ts often difficult, [ff not 
Imposstble, to evaluate exactly. Frequently, one Is not so 
Interested In the exact solution but the behavior of the 
system when one parameter becomes very large or small. This 


is where asymptotic analysis Its needed, 


MACSYMA {fs used to obtain the asymptotic solution of an 
inffinite Integral. Many facllitles provided In such an 
algebralc manfpulatton system can be Illustrated through 
this application. It Is also posstble to show how the 
definite Integration capabllity provided by WANDERER Is 


needed for the successful computation of the results. 
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1. An Asymptotic Analysts Problem and an Outline of the 


Method of Steepest Descent 


Consider the iInfintte Integral 


? ¢ H(t) 
JC) = SE dt, C>0 


where H(t) == -- = Zit 
4 


While tntegrating J exactly may be Imposstble, [ts 
asymptotic behavior as C becomes very large can be 
Investigated. To obtaln the asymptotic expansion of J, the 
method of steepest descent [5,7] will be employed. 
BasIcally, the method of steepest descent consists In 
deforming the contour of Integration In such a way that the 
major contrtbutton to the Integral artses from a small por- 
tion of the new path of Integration. The contribution wil] 
become more and more dominant as the parameter of Interest 


grows. This parameter here ts C. 


The flrst step In thls method ts to find the new path 
of Integration, On a given contour, larger contributlons 


comes from portions where the Integrand ts larger In 
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absolute value and Jess oscillatory. Hence the requirements 
for a desirable path are: (a) the absolute value of the 
Integrand becomes maximum at a pofnt, tl say, on the path. 
(b) The argument of the Integrand ts constant on the path 
near tl. The first requlrement is obvious. The second Is 
essential, for if the phase angle changes even slightly near 
tl, thls change will be magnified by the very large factor C 
resulting in rapid oscillattons of the Integrand and 
therefore negating any possible contribution from the poltnt 
tl. Let UCt1,t2)= RICHCeL4Zt t2)) and VCtl,t2) = ImCHCt1l+31 


t2)). Let p be a potnt where 


dU dU dV dV 


weer Bonner scene ane = 0, 


d tl dq iZ d tl d.t2 


Then p ts certainly a candtdate for tl required In (a) and 
(b). Such a point Is called a saddie potnt. The chotce of 
the name saddle potnt will be made clearer later. The 
Cauchy-Rieman conditions [Imply that criterton for such a 
point Is H'(t)=0. There may be more than one such pofnt in 
the complex t=plane, There are an Infinite number of diffe- 
rent curves which pass through a saddle potnt and satisfy 
(a) and (b). Among them the path along which the integrand 
decrease In size most rapidly ts the best. for the Integral 


J, this means a curve on which U decreases most rapidly. 
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Let H'(b) = 0. It can be shown by use of properties of 
analytic functions that U varles most rapidly on curves 
Im(H(t))=const. (usually called level curves). If H''(b)#0, 
two level curves wlll pass through the polnt b, [tntersecting 


at right angle, as shown tn fig. 1. 


fig. J. 


On one of these two curves, say curve A, RICH(t)) ts minimum 
at b and Increases as t moves along A away from b. On the 
other hand, RICH(t)) Is maximum at t=b on the curve D and 
decreases as t moves away from b on D. Curve A fs called 
the steepest ascent path and D the steepest descent path. 
The point b Is usually referred to as a saddle pofnt. If 
H''(b)#0, b Is a saddle polInt of order 1, If H''(b)=0 and 
H'''(b)#0 of order 2, etc. If the orlginal contour can be 
deformed onto one or a combination of such steepest descent 
paths, then the asymptotic expansion of the given Integral 


can be obtained by a rather rout!ne procedure which 
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Involves: change of varlable of Integration, Inverston of 


truncated power serfes at the saddle potnts and term by term 


Integration. 


Let us outlIne the steps of the solutton procedure as 


follows. 


1) Locate and determine the order of saddle potnts 
of H(T). 
2) Compute UCt1,t2) and V(tl,t2) such that 

H(tl + Z1 t2) = UCt1,t2) + ZI V(tl,t2). 
3) Obtain V(tl,t2) = const. curves which pass 
through the relevant saddle potnts. 
4) ExamIne the V(tl,t2) = const. curves to 
determine whether deformatton of contour can be made 
to curves through the saddle points. 
5) Change the varlable of Integratton. 
6) Express t as a truncated serles In the new 
varfable about each relevant saddle pofnt. 
7) Determine the coefficients In the above series. 
8) Apply Watson's lemma to obtain the first few 
terms of the asymptotic expanston by Integrattng 


term by term, 
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2. Solution Steps 


Presented here are the solution steps of the asymptotic 
expanston of J(C) In the exact sequence as they have heen 
carrfted out using MACSYMA. The lfnes labelled (Cl) are 
Input commands. A command !tne ends with elther a @ sfgn or 
a $ sign. The @ slgn causes results obtafned by executing 
the command ! ne to be displayed In a subsequent lIne 
labelled (DI). A (CI) line together with a corresponding CDi) 
line will be referred to as step I. The $ as an end of 
command I!Ine character suppresses display of results for 
that line. Explanatory texts wlll be Inserted between 
lines. To avotd becoming a user's manual for MACSYMA, 
explanation for the commands used In the solutton will be 
made quite brief. For a more detalled look at MACSYMA the 


reader ts referred to [16]. 
(C1) P:8$ 


P Is a parameter which Is set dependtng on the number 
of terms desired In the asymptotic expanston. By setting P 
to 8, we shall obtaln the first & terms. The reason will 


become evident later. 


(C2) HOT) 23-T##*h/ 4-31 "TQ 


14k 


4 

r 
(D2) H(T):2 - -- = 31 OT 

m 

defintng the functton H(T) 
(C3) DIFFCH(T),T)=0@ 
5 

(D3) -T - 1 = 0 


creating an equatlton H'(T)=0 


(C4) SOLVE(Z,T)@ 


SOLUTION 
SQRT(3) = 3! 
(E4) JT = en nen enna me 
2 
(E5) T= 31 
- SQRT(3) - 3! 
(E6) T = eee e nnn enn nnn 
2 
(D6) (E4,£5,E6) 


The % stgn used In (C4) stands for the last (D3), In 
general a % slgn represents the last expression labelled 
(DI). SOLVE In (C4) fs an Invocatlon of the MACSYMA SOLVE 
program (see [10] for Its capabIlIitftes and lImltattons). 
The roots of H'(T)=0 give three first order saddle potnts,. 
Passing through each of these polnts there wlll be one 
steepest ascent and one steepest descent level curve 


(Im (HCt))= const. curves). The saddle pofnts fn (Eh) and 
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(E6) are the potnts B and A respectively shown In flg. 2. 
The value of H at these three saddle polInts will now be 


computed. 


(C7) HB:RATSIMPCHCPARTCE4, 2)))@ 
- 3 %1 SQRT(3) - 3 

CD73: ee a ae te aims 

What has been done In step 7 Is the computation and 
stmplificatton of H(B). The command RATSIMP causes rational 
stmpliffcatiton which Is essentfally puttIng expresstons to 
be simplified Into the form of one numerator, and one 
denominator and perform all posstble GCD cancellations. The 
command PART allows a user to obtaln subexpressions of an 
expression. PART(E4,2) returns the second part of the 
equation E4 which Is Its right-hand stde. Note that 


‘commands can be nested. 


(C8) HC: RATSIMPC(H(PARTCE5,2)))@ 


(D8) Ps 


(C9) HA:RATSIMPCH(PART(E6,2)))@ 


(D9) ERs aE chr MON 


The potnt T=%1 turns out to be trrelevant because the 
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new path which wlll be determined does not pass through It. 
The next goal Is to obtain curves passing through the saddle 
points (SQRT(3)/2, -1/2) and (-SQRT(3)/2, -1/2) along which 
the Imaginary part of H(T) Is constant and the real part of 
H(T) varles most rapidly, t.e. the steepest paths through 


the saddle poitnts. 


(C10) EXP: EXPAND(H(T1+21*T2))@ 
4 2 2 


T2 3 3T1 T2 
(D10) - --- + 31 TL 12) + eneeeen=- 
4 


3 Tl 
-%1 Tl %T2 + T2 = --- =-%1 Ti 
4 


(C11) V(T1,T2):="COEFF(D10,21)@ 
3 3 
(D11) V(T1,T2):= T1 T2 - T1 %T2 - T1 


By steps 10 and 11, the tmaginary part of H(T1+ 31 T2) 
ts found and given a function name V(T1,T2). The command 
COEFF(Cexp,var,n) computes the coefficlent of var**n In exp. 
In this case, the coefflclient of %1 In (D10) Is exactly the 


Imagtnary part of (D10). 


(C12) UCT1,T2):="COEFFCEXP,%1,0)@ 
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l 2 2 \ 
T2 3 Tl 72 T1 

(D12) WT1,T2):= = --= + --------- + T2 = --- 
m 2 l 


U(T1,T2) Is deftlned to be the real part of 
H(T1 + %1 T2). 


(C13) V(T1,T2)=V(SQRT(3)/2, -1/2)@ 


3 3 3 SQRT(3) 
(D13) TLT2 - Tl T2 - Tl = ----+---- 


Obtained In (D13) Is the equation of a curve, V = 
constant, which passes through the saddle potnt B, 


(SQRT(3)/2, -1/2). This curve wlll be referred to as CR. 


(C14) V(T1,T2}=V(-SQRT(3)/2, -1/2)@ 


3 3 3 SQRT(3) 
(D14) T1 72 = 71 72 -T1 = ----- aan 


(D14) Is the equation of a curve, V=constant, passing 
through the saddle point A, (-SQRT(3)/2, -1/2). Let thls 
curve be CL. The curves GL and CR have to be examined 
carefully by the human problem solver to determing the new 
path of Integration. The manner In which they extend to 
Infinity Is often Important In deforming the contour, 


Asymptotes to these curves can be found easfly. For an 
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ay-=0¥ 
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ie 
B/ElE—-=(2k FON yy BLEIE=(Y¥ PANE D 


2a¥ 


er 


It can be seen tn flg. 2 that by use of Cauchy's 
Integral theorem, the ortg!nal contour can be deformed to 
the contour Ca+Cbh. For the purpose of change of path of 
Integration a new vartable R ts tntroduced which will be the 


parameter of our steepest descent paths. 


(C17) -R**2=H(T)-H(TO)@ 


h 

2 T 1 

(D17) -R = %1T70---- 31T+- 
i" 4 


TO stands for elther of the two saddle pofnts A or B. 
Indeed if TO=A, R ts real If and only If T ts a potnt on Ca, 
for ImCH(T)-HCA))=0 and 0 2 (H(T)<HCA)) only for polnts on 
Ca. The same can be sald about T0=B. T has to be expressed 
In terms of R In order to perform the change of variable. 
Solving for T as a function of R exactly In (D17) fs not 
necessary. What Is needed {ts the first few terms of a power 
serles expansion of T fn terms of R about the polnt T=TO 


(the saddle pofnt). 


(C18)T=T0+DOSUM(JU,1,P, BJU] *Re*Jd)@ 


8 7 6 5 4 
(D18)T = TO +B R +B R +B R +B R +B R 
8 7 6 5 4 
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T Is set to a truncated power serfes In R. The unknown 
coeffictents Bi wlll be solved for by use of (D17). 


(C19) E:SUBSTITUTE(D18, R**2+PART(D17,2))@ 


8 7 6 5 rf 
(p19) - (TO +B R +B R +B R +B R +B R 
8 7 6 5 k 
3 2 4 8 7 
+B R +B R +B R) /4 - 1 (TO +B R +B R 
3 2 l 8 7 
6 5 4 
+B R +B R +B R +B R +B R +B R) 
6 5 4 3 2 


hs 


2 
+ %71 TO +R + 


The expresston (D19) has been given a name E. Before E 
Is expanded In order to collect terms, a few simplification 
rules are defined so as to discard powers of R higher than 8 
In the expanston process. This greatly reduces the 


Intermediate expresston bulge whitch would otherwlse occur, 


(C20) DECLARE (N,NPRED)$ 

(C21) NPRED(X):=IF X > P THEN TRUE ELSE FALSE$ 
(C22) TELLSIMP (R**N,0)$ 

ts THE REPLACEMENT 


A simplificatton rule has been set up so that any R**n 


ts replaced by 0 {tf N > p which ts 8. 
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(C23) TELLSIMPCTO**3,-21)$ 


An addItlonal stmplificatton rule Its deftned so that 


TO**3 shoil be replaced by -%f automatically. 


(C24) RATVARSC(R)$ 
R shall be the main varlable In subsequent rational 


simplificattons. 


(C25) E:PARTCRATCE),1)@ 


(D25) 
8 2 2 
R((C - 12 8B B = 128 B -128B B = 68 ) TO 
1 7 2 6 a) 4 
2 2 
+-C = 1278 B - 24 B B B +( = 24 B B - 12 B ) 8B 
i 6 do 205 1 3 2 4 
2 3 2 2 2 
- 12B B )T0=-4 B B = 12 B B B = 6 8 B 
3 1 5 1 2 4 1 
2 4 
= AZ) BB B -B ) 
1 2 3 2 
7 2 
+R (( - 12 8B B = 12 B B = 128 8B) TO 
1 6 2 5 3. OU4& 
2 2 2 
#00 = 228 B - 24 B B B = 12 B B - 12 8 B ) TO 
1 5 1 2 4 1.3 2 3 
2 3 
- 4B B - 12 B B B-=-4&B B ) 
2 3 1 2 
6 2 2 
+R (€C - 12 B B - 128 B =-=68 ) TO 
1 5 2 4 3 
2 3 
+ ( - 12 B B - 24h B B B =-=4& B ) TO 
a 4 i :2- 3 2 
3 2 2 
- 4 B B - 6 B B ) 
1 3 1 2 


2 
+R (€C - 12 8 B -+ 12 8B 8B) TO 
1 4 2 3 
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2 2 3 
+ ( - 12 8 B - 128 B )TO -4& 8B B ) 
1 2 1 2 
k 2 2 2 4 
+R (C - 12 8B B =- 6B ) TO - 12 8B B TO -8 ) 
1. 3 2 1 2 1 
3 2 3 2 2 2 
+R ( = 12 B B TO - & B TO) +R (4 = 6 B TO ) 
1 ‘2 1 1 


E Is set to the huge expresston above which ts the 
numerator of a truncated expansfon of E3; the denominator Is 


a constant. 


(C26)FOR J:1 STEP 1 UNTIL J > (P-1) 
DO BLOCK CEQ{U]:COEFFC(E,R,J+1)=0,DISPLAYCEQ(uU] ))$ 


This do loop ts used to generated the coefficlents of 
the varfous powers of R In E and store them [In the form of 
equations tn an array RL. Each EQ(JU] will he displayed 


after It {ts set. 


2 

4 - 6 B TO = 0 
1 
2 

- 12 B B TO = & B TO = 0 
1 2 
2 2 2 4 
Ce 12 6B B = 6.3. ) TO = 12-8 B TO - B = 0 


2 
+ ( - 12 8B 
1 


B 


= 12 


(-12B 8B = 128 B 
2 5 


2 3 
- 128 B )T0- 4 B B 
2 3 Z 5 


B B 
2 5 


2 3 
B B ) TO - 4& 8B B 
1 2 1 
2 2 
- 6 B ) TO 
3 
3 
B B =~- & B ) TO 
1 2 3 2 
2 2 
B - 6 B B 3 
3 1 2 
2 
~ 12 B 8B) TO 
3°C«dWd( 
2 
- 12 B B - 12 8 
2 3 
2 8B B B = & B 
1 2 
2 
-~ 12 8B B =68 ) 
3.6CU¢SS 4 
+ ( = 24 B B = 12 8B 
1 3 
2 
- 12 8 BB 
1 2 4& 


Belng consclous of the lack of storage space left for 
subsequent mantpulattons, we get rid of the computation 
history to create some space. 

(C27) KILLCHISTORY)@ 
(C1) SOLVECEQ(1),B8[1])@ 
SQRT(6) 


(El) = meow 
1 3 TO 
SQRT(6)> 
(E2) B Bos aeecen- 
1 3 TO 
(D2) ; (E1,E2) 


The choice of the value for Bl here affects only the 
sense wlth which the steepest decent paths are traversed. 
Thus If we choose (£1), R would vary from INF to MINF on Ca 
and from MINF to INF on Cb. 

(C3) B[1]:PART(E1,2)$ | 


A do loop Is used In (C4) to solve for the remalning 
B's (B2 through B7), 
(C4) FOR J:2 STEP 1 UNTIL J > (P-1) 
DO BLOCK (SOL: EV(SOLVE(EQ(J],B[JU]),EVAL), 
B[U] :PARTCSOL,2),DISPLAY(SOL) )$_ 


The EV command, with EVAL as the third argument, causes 
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The unknown coefficients BI! have been determined, The 
next goal fs to compute 


2 2 dT 
EXP(-C R ) -= dR 
7 dR 


[oe] 


dT 
for TO=A and T0=8,. Odd terms of R In -- do not contribute. 
dR 
Therefore only 4& terms need be [ntegrated. These terms are 
even In R, thus the range of Integrattion can be changed to 


(O INFI. A functlon FNCX) will be defined for carryling out 


this term-by-term integration at a varlable poftnt TO=X. 


(C5) FNCX):=BLOCKC(ANS:C, FOR J:1 STEP 2 UNTIL J > (Pe-1) 
DO ANS: ANS+J*RATSIMPCEV(BlU|,TO=X))* 
DEFINTCZE**(-C#Ra*2) eRee(J-1),R,0,1NF),RETURNCANS) )$ 
(CG) FNC-SQRT(3)/2-%1/2)@ 
1S THE EXPRESSION 
post eles NEGATIVE, OR ZERO 
NEGATIVE@ 
For the purpose of Integration, WANDERER askes the user 
about the sign of €, The contribution from saddle potnt A 


Is given in the next result. 


SQRT(6) SQRTCZPI) 


Coc SOR TCS) 2.3 ie - SORTCE) 


3 SQRTCZPI) (3 SQRTC3) + 23 31) 


4 SQRT(6) (27 31 SQRT(3) + 27) C 
5 SQRT(ZPI) (25 SQRT(3) + 247 31) 


5/2 
5184 SQRT(6) C 
- (35 SQRT(ZPI) (SQRT(G) (672 SQRT(3) + 1816 31) 
3582 SQRT(3) 13077 31 7/2 
w aanneee feemen 3 en eeeene ))/(839808 Cc) 
SQRT(6) SQRT(6) 

(C7) ENCSQRT(3)/2-%1/2)@ 
1S THE EXPRESSION 
-C 
POSITIVE, NEGATIVE, OR ZERO 
NEGATIVE 


SQRT(6) SQRT(ZP1) 
(3 SQRT(3) - 3 31) SQRT(C) 
3 SQRT(ZPI) (3 SQRT(3) - 23 31) 


h SQRT(6) (27 Zt SQRT(3) - 27) C 
5 SQRT(ZPI) (25 SQRT(3) = 247 41) 


5184 SQRT(6) C 


+ (35 SQRTCZPI) (SQRT(6) (672 SQRT(3) - 1816 41) 


3582 SQRT(3) 13077 31 7/2 
bo meneneennnne = neeee-- ))/(839808 ©) 


SQRT(6) SQRT(6) 


The above Is the contribution from saddle potnt B. Now 
the final result Is obtatned by computtIng the followtng 


(-3 J3 %1-3)/8 (3 J3 %1-3)78 
E D6) 


2 (% D7 - 2 %E 


where D6 and D7 are as given above. This expression ts, by 
Inspection, equivalent to 


(-3 V3 21-3)/8 (3 J3 31-3978 __ 
E D D7) 


2 (% 7 + SE 


where D7 Is the complex conjugate of D7. Therefore, 


(C8) 4*RLCZE**HA*D7)@ 
3 SQRT(3) ZI - 3 


(D8) 4 RLCZE a GPa eg NGAP 
(3 SQRT(3) - 3 31) SQRTCC) 


3 SQRT(ZPI) (3 SQRT(3) = 23 31) 


4 SQRT(6) (27 31 SQRT(3) = 27) C. 
5 SQRT(ZP1) (25 SQRT(3) - 247 41) 


5184 SQRT(6) C 


+ (35 SQRTCZPI) (SQRT(6) (672 SQRT(3) - 1816 31) 


3582 SQRT(3) 13077 31 7/2 
to mee neenennnn 2 annenene ))/(839808 C+)? 
SQRT(6) SQRT(6) 
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This ts the first four terms of the desfred asymptotic 
expansion. RL Is not a command of MACSYMA, It {fs used here 


to denote the real part of an expresston. 
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CHAPTER VIII 
CONCLUSIONS AND SUGGESTIONS FOR FURTHER WORK 


The evaluation of defintte Integrals Is a classical 
problem tn mathematics. Great Ingenulty fs frequently 
required, with many Integrals demanding special devices. 

The lack of a suffictently general theory makes evaluation 
of definite Integrals very difficult. {t Is doubtful that a 
theory, comparable In generallty to the Risch Integratton 
algorithm [26] for Indefintte ftntegrals, can be developed fn 
the near future. We have shown that the convergence of a 
class of Integrals of elementary functions {fs recursively 
undecidable. The proof Is In append!x E. The WANDERER pro- 
gram presented here Is a prototype heurfstic computer pro- 
gram for the symbolic evaluation of definite Integrals. It 
contains both general methods such as contour fntegratton, 
residue theory and differentiation with respect to a 
parameter, In addition to quite a few special methods for 
specific types of Integrals. Clues as to which method to 
use for a given problem are obtatned from the Integration 
range and the form of the Integrand. Although many types of 


definite Integrals can be evaluated by WANDERER, It, as 
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almost any other heurlstic computer program, has [ts limita- 
tions. Foremost among these {fs the fact that It Is a pro- 
gram for the evaluation of real defftnite [Integrals of 


elementary functlons. 


It Is hoped that the work reported here may provide a 
starting polnt for new approaches to the evaluation of 
definite [Integrals from the viewpofnt of symbol mantpula~ 
tlon. The advantage of this approach {ts twofold: (1) A com- 
puter can use Integration methods that are too lengthy or 
complicated to be carrfed out by hand; (2) Such a computer 
program contains a collection of powerful methods that can 
Interact with one another and can produce answers to Inte~ 
grals not present In any finite table. For thlis reason, a 
good definite integration program together with other 
facilities provided [tn an algebraic mantpulation system can 


be very useful to applied mathematIcltans. 


DELIMITER ts a rather sophisticated program for comput- 
Ing limits of elementary functions. Such a program has been 
shown to be a very useful tool [In an algebraic manIipulation 
system. The method of comparing orders of Infinity used In 
DELIMITER Is an Important concept which Is useful [fn places 


other than the computation of limits. 
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A bastc assumption of DELIMITER ts that the operator 
LIMIT ts distributive over the operators +, -, *, and **, 
I.e., rules such as the limit of a sum Is the sum of the 
limtts hold. These rules are valid only If the subproblems 
thus generated produce answers which do not lead to an 
Indeterminate form. Of course there Is L'Hospital's rule 
which can be applled to some Indeterminate forms, hut 
generally only those Involving only the operator * or #*, 
Although many aspects of the Indeterminate form CINF-INF) 
have been cons!dered. In chapter 2, the algorithm for tts 
determination Is not complete. A powerful serfes expansion 
program would be helpful In some cases. Yet, such a pro- 
gram can not solve all problems. Consider for Instance, the 
limit problem 

2 2 
LIMITCSIN X + COS X,X, INF). 

To obtaln such limfts, an algebrale mantpulatton system 
must be able to detect all constant Identities, This Is not 
possible for the set of all expressfons [22]. It may be 
possible for a proper subset of all expressions. For 
example, many trigonometrical [dentitles dlsappear If all 
trigonometric functlons are radically transformed [tnto sums 


of complex exponentials. This Is a basic problem of great 


165 


practical and theoretical [Interest [6]. Such sImplification 


capabIlittes would be of use to many other programs as well. 


DELIMITER can be augmented by Increastng the types of 
functftons tt can handle which may Include functions defined 
by f[ntegrals. Another possible area of research fs the 


automatic determination of superfor and tnfertor limits. 


The fact that WANDERER can evaluate many non-trivial 
definite tntegrals does not mean that It can compete with an 
expert human [Integrator yet. For one thing, a mathematictan 
can usually construct a function of a complex varfable and a 
suitable contour for evaluating different Integrals. 
WANDERER cannot form a contour based on analysts of a given 
real Integral. It simply selects from the cases known to 
Tt. WANDERER would be much more powerful {ff [ft could 
determine, for a given Integral, whether the method of con- 
tour integration and res!tdue theory were applIcable, and, If 
It were, evaluate the Integral by forming a suftable I[nte- 


gral around a closed contour, 


The evaluation of contour fntegrals by residue theory 
usually requlres the solutton of algebrafc or transcendental 
equations. WANDERER uses the SOLVE program of MACSYMA for 


this purpose. SOLVE has Its lImitatltons and fts Improvement 
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is also a posstIble area for further work. 


Suppose one wants to evaluate a closed contour Integral 
by restdus theory. An Interesting general question Is: what 
knowledge about the Integrand or the functfons used tn form- 
Ing it Is necessary. We think the follow!ng are essential: 
(1) evaluation; 

(2) differentiabillty and derlvatives; 
(3) singularittes; 


(4) asymptotic behaviour, 


A natural extenston of the work here Is Integration 
over arbitrary user-~specifled contours. This should not be 
difficult to do, except for the lack of notation. 
Specifically, we must specIfy an arbitrary contour to a com- 
puter and devise a general data structure for use In 
representing contours. For ftnstance, consider the 


specification of the following Indented contour, 


a XL 
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One obvious way Is to specify a contour plecewlse. 
Each piece of curve has a parametric form and a starting and 
ending value for the parameter. There may be other 


approaches. 


Further work In thls area of symbol manipulation might 
Include: 
(1) Destgn of computer algor!ithms for testing convergence 
and divergence of Integrals, 
(2) Summation of Infinite serfles by restdue computations, 
(3) Investlgation of algorithms for definite. Integratton of 
_spectal functtons. 
(4) Evaluation of multiple definite Integrals. 
(5) Application of defftnite Integration programs [fn solution 


of differential equattons. 
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APPENDIX A 
A FLOWCHART OF DELIMITER 


A flowchart {fs presented In the next page which detalls 
the flow of control of DELIMITER, Listings of the programs, 
written In the LISP programming language, may be obtalned 


from the author. 


The routine LIM, appearing In the flowchart, fs a 
program which appltes the 'trivlal' limlt rules, makes use 
of subroutines to compute Ifmits of sums, products, powers, 
and the functions SIN, COS, TAN, LOG, SINH, COSH, TANH. It 
calls LIMIT recursively and makes use of L'Hospital's rule 


and other routines when needed, 
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FE __WRETURN ANSWER] 
LIMITCE(X),X,L) L'HOSPITAL'S RULE RETURN ANSWER 
FAILED? 

ETARD : 
E, X, 
OF X? 
= 
R E 


COMPARE ORDERS 


es 
LYEHOSPITAL'S RULE N AND D 


FE INDEPENDENT 


N AND D CONTAIN 
FUNCTIONS IN THE 
FORM e**R(X) 

WHICH TEND TO INF? 


ETURN 


E RATIONALL F 
IN xX ? 
T 
RETURN LIMCNL/DL) 
RATIONAL NLAND DL= 1 
CUNCTION ND D NFINITY?| 
ROUTINE 


NL=INF OR MINF 
AND ? 
DL=INF OR MINE 


N-NUMERATOR 
OF E 
D-DENOMINATOR 


RETURN _IN 


pL=0?} --@ 
‘ RETURN MINF 
RETURN F 
LIM 
F on NL=0 AND DL=07 
N_AND DARE] ¢ RETURN INFINITY 
RP-FORMS? 
a 


NL=LIMITCN, X,L 
RP=FORM ROUTINE]  (OLeLIMIT(D,X,L) 
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APPENDIX 8B 
A PROOF 


Infinite [tntegrals of 
n 
lk Z om 
f(Z) = ZE Z 
where n > 0, k # 0 real, R1I(m) > =-1, and n = RI(m) > 1, 


have been discused [tn chapter 4-3.2. Some results dertved 


there depends on the proof of (1) and (2) below. 
The objective here ts to supply the proof of 


1) inte [ F(Z) dZ = 0, 
é¢-—" OF Ce 


and 


2) diode f f(Z) dZ = 0. 
R*+0 JoR 


where CR and C are the clrcular contours gtlven below. 


y 


Ce x 


CR 


Let 
te 
Z=re and m= a+ tb, 
then 
fem dZ = 
c 
m+1 man n Ite 
Ir EXP(T kre +i(m+t1)8@) de 
ro) 
atl wan n 
(i) <r EXP( = k SIN(n @) r - bo) deo =M 
- wv 


For 0 se S %PI/(2 n) we have 


n 

EXP( - k SIN(n eo) r ) = 1. 
Thus 
Since a + 1 > 0, this completes the proof of (1). 


Now from equation (1) we have 


%/2. 5 iia 
EXP( = k SIN€e) r = ——) do 
°o 


n 


atl 


n 


eles: 


ee a n ZPt b 

—~m(] = EXPC = k fr = 

n 2n 
Zkonr + SPI b 


in 


)) 


The fact that m-a-1 > 0 completes the proof of (2), 
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APPENDIX C 


PERFORMANCE OF WANDERER 


AND DELIMITER 


As a measure of the performance of WANDERER and the 
limit programs, selected problems have been timed. The time 
sharing system used ts the ITS of the Artificial 
Intelligence Laboratory at MIT which uses a PDP-10 computer 
with a memory cycle time of about 2.75 microseconds. Time 
used for parstng the fnput string and display of the 
computed result has been excluded fn order to obtaln an 
approximation to the time actually spent Inside WANDERER or 
DELIMITER. Garbage collection fn the LISP [17] system tn 
which MACSYMA fs written may take place during a 
computation. Although garbage collection Is a slow 
procedure, {ft Is only falr to regard [ft as part of the 
computation process belng carrled out. Therefore, the time 
It requtred has been [Included tn the timing experiments, 

The results are put In the form of two tables. An * 1s used 
in the tables to Indicate computation requiring LISP garbage 


collectton. 


17h 


TABLE 3 


PERFORMANCE OF OELNIETER 


LTALt PRriatees seared Result Time tn 
equivalent to Sec, 


LIMIT x 


A+ O+ 


aed 
[any 
Ce 
haw 


LIMIT (lei/x2_ Ee ee) 


A7% +00 


LIMIT eee INE Louis 


Kine +00 


LEM, “Sees aes oe co a an on INF tS 


pee SADT! 
ron ae ee a 


tenets od MEME Bey 


~ 
hm) 
‘4 


7 Beas “ 


TABLE II 


PERFORMANCE OF WANDERER 


Integral computed 


Result 
equivalent to 


0 Cprinctpal) 


LOG(2) 


3 A LOG(3) + (B +.3)%PI 


20 SINCZPI/20) 
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Time fn 
Sec, 


2.02 


7.42 


11,32 


1.08 


3 D -5 SQRT(X) D 

X CE dX C2016 FE Df P8125 
oO 

ao 2 

COS (X) 

“sera aX 2P1/28 
lo 2 

xX + 49 

e dX 

SSeS SSeS ai ®PI/SQRT(2) 

SQRT(X) (2 X + 1) 
fe] 

~o 1/3 

x aX 


we nn ene DIVERGENT 


°C SORTOCX) * 3) 
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0.95 


6.62 
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APPENDIX E 


AN UNDECIDABILITY RESULT 


Let Sl = {P(X1,X2,...,Xn)} be a set of polynomials with 
integral coefficients In Xl, ... , Xn and S2 a set of 


functions F of the form 


F(X1, aon ¢@ Xn) 
2 2 ” 2 2 
= (ntl) (P (X1, we. , Xn) +). (SEN MXJIKJ (XL, 22. 4 Xn))-1, 
Jel 
where Kj Is the domInating functfon |6| for 
d 2 


— Pp, 


dX j 
Richardson has shown |{22,6] that 


Lemma 1. For F in S2 the predicate "there exists an n-tuple 
B of real numbers such that F(B) < 0" ts recurstvely 


undecfdable,. 


Lemma 2. If FCB) < 0 for some ntuple B real numbers then 
there exists an n-tuple A of nonnegative [ntegers such that 


P(A)=0 and therefore FCA) = <1. 


Lemma 3, Let 


R(X) = 


Then for any real numbers Al, 


there exists b > 0 such that 


thCbd-ALI<E, [h€glb))-A2|<€, 


By use of these Lemmas we 


Corollary. For any G ina set 


GUA SBCA Dy ROR I Dy ware gy 


the predicate 


0" is recursively undecidable, 


Proof: Suppose the predicate 


we have 
1) 


real numbers Al, ... , An such 


2) If there exist real numbers 


FCAl, «2.3 , An) < 0 then (by 


nonnegative integers Bl, ... , 


FCB1l, ... , Bnd=-1,. Thus, 


that G(c) < 0. Since 


G(X) Is continuous, this 


X SINCX) and g(X) = 


zoe , 


"there exists a real 


Is recurstvely decidable, 


there exists a real 
G(X) can be 


Implles that there exist at 


3 
X SINCX ) 


An and any 0 ¢€ ¢ |] 


oe 6 


Phe ss. .80b) 53399 An Ces 


can show 


of functtons of the form 
hCe(...€20X))...))) + 1/2, 


number t such that G(t) = 


then 


If there exists a t such that G(t)=0 then there exist 


that FCAl, ... , An) ¢ 0; 


Al, ... , An such that 


lemma 2) there exist 


Bn such that 
number c such 
Jarge and positive and 


such 
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that G(t)=0. 


This contradicts lemma 1. 


From the definition of the Functlon F, one can see G is 
always -1/2 and G {ts large and pos!ttive except at the 


vicinity of a finlte number of pofnts where G {ts negative. 


Theorem: The convergence of a set of ftntegrals of the form 


@ dX 
[ 2 2 
—20(X +1)G (X) 


Is recursively undecidable. 


Proof: This tntegral Is convergent If and only tf G(X) has 


no real zero, 
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